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Abstract 

We summarize how future measurements of electromagnetic (EM) probes at the 
Relativistic Heavy Ion CoUider (RHIC), in connection with theoretical analysis, 
can advance our understanding of strongly interacting matter at high energy den- 
sities and temperatures. After a brief survey of the important role that EM probes 
data have played at the Super Proton Synchrotron (SPS, CERN) and RHIC to date, 
we identify key physics objectives and observables that remain to be addressed to 
characterize the (strongly interacting) Quark-Gluon Plasma (sQGP) and associ- 
ated transition properties at RHIC. These include medium modifications of vector 
mesons via low-mass dileptons, a temperature measurement of the hot phases via 
continuum radiation, as well as 7-7 correlations to characterize early source sizes. 
We outline strategies to establish microscopic matter and transition properties such 
as the number of degrees of freedom in the sQGP, the origin of hadron masses and 
manifestations of chiral symmetry restoration, which will require accompanying but 
rather well-defined advances in theory. Increased experimental precision, an order of 
magnitude higher statistics than currently achievable, as well as a detailed scan of 
colliding species and energies are mandatory to discriminate between theoretical in- 
terpretations. This increased precision can be achieved through hardware upgrades 
to the large RHIC detectors (PHENIX and STAR) along with at least a factor of 
ten increase in luminosity over the next few years, as envisioned for RHIC-II. 



Key words: direct photon, isolated photon, quark-gluon plasma, relativistic heavy 
ion collision 

PACS: 25.75.-q, 25.75.Nq, 12.38.Mh, 13.85.Qk 



Preprint submitted to Elsevier Science 



25 April 2008 



Contents 

1 INTRODUCTION i 

1.1 Toward Discovery and Characterization of Hot and Dense QCD 
Matter i 

1.2 Where EM Probes Are Unique or Very Important 

1.3 Theory vs Experiment at the CERN-SPS i 

2 THEORETICAL CONSIDERATIONS AND PREDICTIONS M 

2.1 Objectives and Framework [l^ 

2.2 QCD Lattice Results [13 

2.3 Sum Rules [l4 

2.4 Temperature and Degrees of Freedom [l3 

2.5 In-Medium Spectral Functions below and above Tc [13 

2.6 Low-Mass Dilepton Spectra [13 

2.7 Chiral Symmetry Restoration [33 

2.8 Electromagnetic Signatures of the Color Glass Condensate [13 

3 OBSERVABLES: STATE OF THE ART M 

3.1 Low-Mass Dileptons [13 

3.2 The 0(1020): Hadronic vs Leptonic Decays and V2 [41 

3.3 Direct Photon Spectra [13 

3.4 Direct Photon Raa \M 

3.5 Photon Azimuthal Asymmetries (Elliptic Flow) [47 

3.6 Electron Raa and Flow [47 

3.7 Direct Thermal Photons via Low-Mass Dielectrons [13 

4 THE FUTURE M 

4.1 Detector upgrades [53 

4.2 New Measurements \5U 



2 



4.3 High Statistics and Energy/Species Scan 

4.4 From Observables to Pertinent Theoretical Issues 
5 SUMMARY AND CONCLUDING REMARKS 
References 



1 INTRODUCTION 



1.1 Toward Discovery and Characterization of Hot and Dense QCD Matter 



One of the key goals of ultrarelativistic heavy-ion colhsions is the creation of 
hot and dense strongly interacting matter that: (i) resembles the conditions in 
the early universe; (ii) can be related to the phase diagram of the underlying 
theory (Quantum Chromodynamics (QCD)); and (iii) enables the discovery of 
new phases. Recent surveys of results from the first three years of data taking 
by the four experiments at the Relativistic Heavy Ion Collider (RHIC) [1,2,3,4] 
are unanimous in their conclusion that a new form of matter has been created. 
This matter is very dense, opaque and exhibits a high level of collectivity 
which has largely been attributed to the expansion of a partonic phase. It is 
inconsistent with naive expectations based on a weakly-interacting (gas-like) 
Quark-Gluon Plasma (wQGP), while it is best described in terms of a so- 
called strongly interacting QGP (sQGP) constituting an almost perfect fluid. 
Thus, a discovery has been made in a qualitative sense, but the properties of 
this new state of matter remain under intense debate. 

A closer look indeed reveals that we are still quite far from a coherent and 
quantitative description of the sQGP at RHIC. The wealth and precision of 
new data from Run-4 and Run-5 [5] is a first step in this direction. Measure- 
ments of previously inaccessible signals, such as semileptonic electron-decay 
spectra, J /ijj production and three-hadron correlations, and improvements on 
the statistical and systematic errors, as well as the range of the existing data, 
have been achieved. In addition, new and ongoing analyses of SPS data, most 
notably from NA60 [6,7] (including low- and intermediate-mass dileptons and 
J/ip production), have reached unprecedented levels of precision that now 
can distinguish between model predictions which were consistent with earlier 
data sets. Furthermore, issues have been raised again, that seemed to have 
essentially been settled a few years ago, e.g. the energy-loss mechanism of 
jet quenching. It therefore appears fair to say that whereas the existence of 
the new form of matter has been established, we neither understand its micro- 
scopic properties nor deduced from it convincing signals of the phase transition 
itself. While theory has also made substantial progress in the last few years 
by moving from more signal-specific explanations to a coherent description by 
connecting different phenomena and improving theoretical tools, a widely ac- 
cepted "grand scheme" encompassing both bulk and microscopic components 
has not yet been realized. As ideas get refined (largely steered by data), dif- 
ferent theoretical predictions often approach each other, thus increasing the 
demand for higher quality measurements to differentiate between them. 
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Table 1 

Matrix of detector or accelerator upgrades vs physics measurements. X = upgrade 
critical for measurement. O = upgrade important for measurement. 

After very successful operation of RHIC0 and its detectors over the first five 
years of data taking, we are approaching a point where further progress re- 
quires improved experimental capabilities. Table 1 lists important physics top- 
ics that are either beyond our current reach or would be significantly enhanced 
with the indicated detector and luminosity upgrades. These should be put into 
context with overarching questions in the investigation of the sQGP and chi- 
ral/deconfinement transitions, emphasizing those for which electromagnetic 
probes are particularly relevant, including: 

• What are the temperatures and corresponding system sizes of thermalized 
matter at its early stages? 

• How do hadron properties change in hot and dense matter and how are 

^ By 2004 (Run-4) the accelerator exceeded design luminosities by a factor of 2.5 
both in Au+Au and pp collisions. RHIC collided 4 different systems including asym- 
metric d-f-Au interactions, collisions at 6 different center-of-mass energies, among 
others at 19 GeV (Au+Au) and 22 GeV (Cu+Cu), establishing the first overlap 
with the CERN SPS. In 2007 (Run-7) it delivered 3260 ^b"^ integrated Au+Au 
luminosity at 200 GeV, in 2008 (Run-5) the first tests with 9 GeV and 5 GeV per 
nucleon (below injection energy) are planned with Au+Au. 
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Fig. 1. Left panel: direct photon excess ratios measured by PHENIX in ^/s^ = 200 
GeV Au+Au collisions [8]. The curves are the ratios predicted by NLO pQCD [9]. 
Right panel: integrated Raa for photons and vr^s as a function of the number of par- 
ticipants, A'part) measured by PHENIX in ^/s^ = 200 GeV Au+Au collisions [8]. 



hadron masses generated? How do the medium modifications depend on 
temperature and net/total baryon density? 

• Can we deduce signatures of Chiral Symmetry Restoration (xSR) and, if 
so, how is it realized? 

• What are the relevant degrees of freedom in the sQGP? Does it harbor 
bound states and/or resonances? 

• How do the medium properties change if the net baryon density is increased? 
Can a QCD critical point be found? 

• How does the system reach (local) equilibrium on the apparently short time 
scales required by hydrodynamics? 

• What exactly causes jet quenching? 



In the remainder of this introductory section, we will first elaborate in some- 
what more detail on the features and strongholds of EM probes in heavy-ion 
collisions (Sec. 1.2). In particular, we will discuss the achievements and short- 
comings of the EM probes program at the SPS (Sec. 1.3) which will help 
us to sharpen the case for future RHIC measurements. In Sec. 2 we lay out 
the theoretical framework for describing and interpreting EM probes. We will 
also provide predictions for observables and formulate strategies for deducing 
information on key properties of the medium that may not follow from an im- 
mediate interpretation of experimental data. In Sec. 3 we give a brief overview 
on the current status of EM observables at RHIC. In Sec. 4 we then focus on 
how future detector and accelerator upgrades can be used and geared toward 
answering the above questions. In Sec. 5 we reiterate the main points of this 
document. 
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1.2 Where EM Probes Are Unique or Very Important 

Electromagnetic probes are real, 7, and virtual photons (dileptons), 7* l^l^ 
where / = e and /x. They are emitted from the entire reaction volume through- 
out the evolution of a heavy-ion collision, from first impact in primordial 
(hard) NN collisions until decays of long-lived hadrons long after strong inter- 
actions have ceased. Once created, most of them leave the interaction volume 
unchanged due to their negligible final-state interaction with the strongly- 
interacting medium. Following the usual terminology we define direct pho- 
tons as those that are not decay products of final state hadrons. Accord- 
ingly, direct photon spectra are defined as the spectra remaining after subtrac- 
tion/eliminatioiH] of photons coming from final-state decays ("background" 
or "decay photons"), whereas in measured dilepton spectra pairs from final 
state decays are often included and their estimated contribution shown sep- 
arately ("hadronic decay cocktail"). The broad category of "direct" photons 
is then often subdivided according to their source such as "prompt" photons 
(from hard scattering) "thermal" photons, "jet conversion" photons, etc. 

Unfortunately, the same property that allows photons to escape, a « a^, 
also leads to major experimental challenges: low rates and large backgrounds 
from the above-mentioned late hadron decays {e.g. Dalitz-decays for dileptons 
and 7r° — i> 77, 77 — > 77, etc. for real photons). 

The fact that direct photons from initial hard scatterings escape the system [8] 
is demonstrated in Fig. 1 where the measured direct photon excess raticH] at 
sufficiently high transverse momentum, p-r, is shown to be consistent with 
next-to-leading order (NLO) perturbative QCD (pQCD) calculations, even in 
the most central Au+Au collisions. 

The interesting physics we are primarily concerned with here is at significantly 
lower pt and dilepton invariant masses, M than inferred from Fig. 1. Thus, in 
this review, we stress the importance of precision data at masses and momenta 
below ~ 3 - 4 GeV. 

Electromagnetic probes are unique. 

• They give direct access to the in-medium modifications of hadronic states 
(the p(770), u;(782), and 0(1020) vector mesons) via dilepton invariant-mass 
spectra, which can illuminate the nature of hadron mass generation and thus 
the origin of ~ 98% of the visible mass in the universe as well as related 
changes in the structure of the QCD vacuum, including chiral symmetry 

^ This can be done in a number of ways like statistical subtraction, "tagging" , etc. 
^ The ratio of inclusive photons over hadron decay photons. For more details on 
the direct photon excess ratios and Raa see Sec 3.4. 
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restoration. Detecting novel nonperturbative (resonant) correlations in the 
sQGP at masses above ~ 1 GeV is also possible. 
• They can be used to infer the temperature of the system during its hottest 
phases via direct thermal photon and dilepton radiation. In addition, HBT 
interferometry of thermal photons offers the cleanest measurement of early 
system sizes. 

Indirect consequences of the above studies include insights into the mecha- 
nism of rapid thermalization, e.g. via resonance formation in the sQGP. The 
effective degrees of freedom can be determined if a temperature measurement 
is complemented with independent information on the energy or entropy den- 
sity. Although not unique, EM probes are valuable for (i) disentangling the 
energy loss mechanism of jets (jet quenching) by inferring photon radiation off 
energetic quarks and establishing the jet energy scale in 7-hadron jets and (ii) 
obtaining complementary information on (early) matter flow by investigating 
photon elliptic flow, f 2, which, in turn, discriminates between thermal photons 
and those radiated off jets. 

In the following sections, we will elaborate on these statements together with 
the requirements to measure associated observables with sufficient accuracy. 



1.3 Theory vs Experiment at the CERN-SPS 



Before beginning a more detailed discussion of the theoretical aspects un- 
derlying a comprehensive electromagnetic probes program for future RHIC 
runs, we briefly review the main achievements at the SPS prior to the new 
NA60 data [6,7], discussed later. An example of one theoretical approach to 
electromagnetic emission simultaneously applied to semi-/central Pb-|-Au and 
Pb-|-Pb data is compiled in Figs. 2 and 3. The data suggest a common thermal 
source with an initial temperature of Tq ~ 210±30 MeV and a lifetime of about 
12 ± 3 fm/c with a thermal freezeout temperature of around 100 — 120 McV. 
Importantly, the low-mass dilepton enhancement in the CERES e^e^ data [10] 
(upper panel of Fig. 2) requires substantial medium effects over a sufficiently 
long lifetime, on the p spectral function. However, decisive discrimination be- 
tween a dropping-mass scenario [11] and a substantially broadened spectral 
function [12] was not possible since the QGP contribution is small, around 
10 — 15%. The sensitivity to QGP radiation increases in the intermediate- 
mass dimuon spectra of NA50 [13] (lower panel) where the observed factor 
of two excess over the baseline charm and Drell-Yan sources can be reason- 
ably accommodated with thermal radiation containing a 30 — 50% [14,15,16] 
QGP component, the main evidence for the Tq quoted above. A very similar 
decomposition is found in the WA98 direct photon spectra [17] at ~ 2 GeV 
(upper panel of Fig. 3) where q is the real or virtual photon momentum. The 
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tra [17,18] compared to thermal emis- 
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recently published low-momentum data [18] (lower panel), extracted via pho- 
ton HBT methods, are not easily reconciled with theory [19,20]. The inclusion 
of soft Bremsstrahlung off vrvr and nK scattering appears to improve the sit- 
uation [21]. The low-momentum yield is essentially emanating from the later 
stages of the fireball evolution and thus it is proportional to the total fireball 
lifetime. 



The SPS electromagnetic probe program was incomplete for a number of rea- 
sons ranging from insufficient statistics to draw conclusions to the lack of 
experimental cross checks. A partial list of open questions is given below. 

(a) There was no decisive discrimination of in-medium p modifications (al- 
though the recent NA60 data [6,7] has improved the situation, albeit in a 
smaller system, see below). In particular, no systematic excitation function 
was obtained through an energy scan. The only CERES/NA45 low-energy 
run at ii^iab = 40 AGeV [22] indicated an increase over the enhancement at 
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Fig. 6. Temperature 
dependence of the 
quark chiral conden- 
sate, {qq){T), computed 
in unquenched lattice 
QCD [41]. 



160 AGeV but with large errors. 

(b) The Cronin enhancement in the primordial pQCD photons, essential for an 
accurate assessment of the thermal photon yield in the qt — '2 GeV region, 
was not determined. 

(c) Charm dileptons were not explicitly identified, hampering the temperature 
extraction in the intermediate-mass dimuon spectra. However, NA60 made 
substantial progress in an intermediate-size system [23]. 

(d) There was no experimental redundancy to cross-check measurements and 
create competition. 

(e) The QGP contribution to EM probes at the SPS is perhaps too modest in 
principle to reveal itself in a significant way. 

(f) The level of early thermalization at the SPS may be insufficient to justify 
thermal approaches in the early phases, as may be indicated by the lack of 
elliptic flow in hadron data relative to hydrodynamic models above pt ~ 
IGeV. 

We will argue below that all of these issues can be (are) overcome at RHIC-II. 



2 THEORETICAL CONSIDERATIONS AND PREDICTIONS 

2. 1 Objectives and Framework 

Ample evidence for (early) thermalization of the matter formed in (semi- 
/ central) y/s^ = 200 GeV Au+Au collisions from current RHIC data pro- 
vides the necessary prerequisite to validate the study of the QCD phase dia- 
gram. So far, most of the deduced features pertain to bulk matter properties, 
including large energy densities, e, well above the critical density extracted 
from lattice QCD [24]. 
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In the following subsections, we will sketch a theoretical framework that pro- 
vides the basis for gaining decisive new insights in three areas: temperature 
measurements; vector meson spectral functions; and the origin of mass and 
chiral symmetry restoration. 

Explicit temperature extraction has so far been restricted to hadro-chemical 
and -thermal freezeout, see Fig. 4. The goal is to estabhsh early temperatures 
well above Tc, feasible with photons and dileptons at a typical energy/momen- 
tum/invariant-mass scale of 1 < M < 3 GeV. Together with information on 
the energy and/or entropy density, the effective number of degrees of freedom, 
dof, may be accessed since e = l^dofT"^ and s — ^dofT^ [26,27,28]. 

Dilepton invariant-mass spectra are invaluable and unique means of extracting 
microscopic information on the constituents of the medium. Modifications of 
low-mass vector mesons {V = p, uj, and 0) in hot and/or dense hadronic 
matter have been extensively studied theoretically (see Refs. [29,30,31,32,33] 
for reviews), largely triggered by the intriguing excess radiation observed at 
the CERN-SPS [7,10]. At RHIC, for the first time, these measurements will 
be performed in an environment that is close to net baryon-free, which will 
provide important tests of the relevant mechanisms underlying the predicted 
medium effects. In addition, the higher initial temperatures achieved at RHIC 
will more directly access radiation from the (s)QGP and thus be more sensitive 
to suggested resonance/bound-state formation above Tc in the vector channel 
[34,35,36,37] (see Fig. 5), relative to more indirect probes of resonances, e.g. 
quantitative analyses of energy loss [38] and hadron elliptic flow, especially in 
the charm sector [39]), or systematics of charmonium regeneration [40]. 

A more ambitious goal is to infer signatures of the phase transition, requiring 
the study of order parameters. For the chiral transition, these are e.g. quark 
condensates, see Fig. 6, the constituent quark mass and the pion decay con- 
stant. While none of these are (viable) observables, the condensate structure of 
the in-medium ground state is encoded in its hadronic excitations. It is thus 
necessary to establish connections between the in-medium vector correlator 
measured in dilepton spectra and order parameters. 

All of the above are inevitable consequences of QGP formation, albeit mostly 
nonperturbative in nature. While data interpretation will require the applica- 
tion of phenomenological approaches, well-defined links to finite-temperature 
lattice QCD computations and symmetry constraints are essential for the de- 
duction of meaningful results. 

We emphasize that thermal production rates for photon and dilepton spectra 
can be cast into a uniform theoretical framework according to 
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where the key quantity is the retarded electromagnetic correlation function, 
Hem [42,43]. In the vacuum, this function can be measured in e~^e~ annihila- 
tion and decomposes into two regimes: at masses above M = ~ 1.5 GeV, 
the strength of the EM spectral function, Imllg^, is rather accurately deter- 
mined by perturbation theory, i.e. annihilation into qq pairs with little impact 
from subsequent hadronization. At low mass the cross section is saturated by 
the light vector mesons p, cu and 4> (vector dominance), i.e. nonperturbative 
resonance formation: 



imn™^ 

em 



(M) 



E 

v =p,ui4 
-Ml 1 

127r ' 



9V 



ImDviM) 



1 + ^^4r^ + . . .) A^c E e2 M > 1.5 GeV . 



M < 1.5 GeV, 



(3) 
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For thermal dilcpton emission this implies that, on one hand, the low-mass re- 
gion carries information on dynamical medium effects with a relative strength 
of the vector mesons approximately given by 11:1:2 for p : uj : (f), reflect- 
ing the vector-dominance couplings, rriy/gy, or, equivalently, the dilepton 
decay widths. On the other hand, at intermediate mass a reasonably con- 
trolled emission strength provides the basis for probing the temperature. The 
temperature and volume dependence in space-time integrated dilepton spec- 
tra combine such that the prevalent low-mass contribution originates from 
temperatures around and below Tc [44]. At higher masses and energies, the 
exponential sensitivity of the Bose factor strongly biases contributions toward 
high temperatures. Similar considerations apply to thermal photon transverse- 
momentum spectra, corroborating the feasibility of realizing the three basic 
objectives listed above. Note that the leading-order contribution to the dilep- 
ton rate at sufficiently large invariant mass, M > 1.5 GeV, is C(q;°) while a 
nonzero photon rate requires processes of at least 0{as). 

In rough accord with our physics objectives, we adopt the following classifica- 
tion of regimes in dilepton invariant mass, Mu, or photon transverse momen- 
tum, qt: 

• the low- mass region (LMR), Mu < 1.1 GeV (vector meson decays); 

• the intermediate-mass region (IMR), 1.1 < Mu < 3 GeV (continuum radi- 
ation, QGP emission, resonances?); 

• the high- mass region (HMR), Mu > 3 GeV (primordial emission and heavy 
quarkonia) . 
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Fig. 7. Thermal dilepton production rates in the QGP as evaluated in quenched 
lattice gauge theory at zero three-momentum (dotted and dash-dotted line) [45] . The 
solid and dashed lines are pQCD calculations of the rate in the Born approximation, 
qq annihilation to leading order, 0{a.^), and with hard-thermal-loop corrections [46] 
respectively. 

2.2 QCD Lattice Results 

Thermal dilepton rates have been studied in lattice QCD within the quenched 
approximation [45]. The computation of the finite-temperature Euclidean cor- 
relators in the vector channel is supplemented by a transformation into the 
time-like regime using the maximum entropy method, after which the dilepton 
rate follows from Eq. (2). The results for zero three- momentum at two differ- 
ent temperatures above are compared to perturbative calculations [46] in 
Fig. 7. The lattice rates are quite comparable to the perturbative ones at high 
energy, while the enhancement at u; ~ 5T could be related to resonance for- 
mation in the QGP. Toward small energies, the IQCD rates drop substantially 
and deviate markedly from the perturbative calculations, which, in the HTL 
approximation, even diverge. Since the HTL calculation is of nontrivial or- 
der in as, the HTL divergence is closely related to a non- vanishing thermal 
photon rate. Recent lattice computations of the vector correlator at finite 3- 
momentum [47,48] confirm that the rates in Fig. 7 lead to vanishing photon 
emission. If this seemingly unrealistic feature is induced by lattice artifacts, the 
decrease of the dilepton rate toward small cu may not hold, leading to better 
agreement with HTL extrapolations, consequently affecting the interpretation 
of the "resonance" structure. The ultra-soft limit of the EM correlator can be 
further studied by its relation to the electric conductivity, cTem, via the Kubo 
formula, 

(JemiT) = y ^ImnJ^(go, q = 0; //b, T) . (4) 

The correlators found in Ref. [45] correspond to a vanishing conductivity. A 
different method employed in Ref. [49] leads to a finite, large value of cXem, 
o'em{T) ~ 7J2^lT, see Ref. [48]. It would be very interesting to compare this 
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result to the conductivity underlying the soft photon emission rates used in 
the description of the low-momentum WA98 data [18], see the lower panel of 
Fig. 3. Along these lines, (Xem has been evaluated in a low-temperature pion 
gas within chiral perturbation theory 



Another quantity of interest is the susceptibihty associated with conserved 
quantum numbers, i.e., derivatives of the thermodynamic free energy with 
respect to a chemical potential. 



XX = -^-r • (5) 



These susceptibihties have been evaluated in lattice QCD for quark and isospin 
chemical potentials [51]. When extrapolated into the finite- /Xg plane, the isoscalar 
quark susceptibility develops a maximum while the isopin susceptibility re- 
mains monotonic. Since the susceptibilities can be related to the space-like 
static limit of the corresponding correlation function in the cu and p channels, 
valuable information on the soft part of the spectral functions may be in- 
ferred [52] or at least tested for a given model. Whether the lattice results for 
the isoscalar susceptibility are related to an in- medium reduced u mass e.g. 
via a (J — a; mixing mechanism close to the QCD critical endpoint, remains to 
be seen. Interestingly, recent experiments on nuclear cu photoproduction have 
provided evidence of cu mass reduction in nuclear matter [53]. A similar line 
of reasoning for the p suggests that its mass is not much affected approaching 
the critical point [54] . 



2.3 Sum Rules 



Besides direct IQCD calculations of dilepton and photon rates, unanticipated 
in the immediate future, model-independent information and effective-model 
constraints are encoded in (energy-weighted) sum rules relating (integrated 
moments of) spectral functions to vacuum expectation values of composite 
quark and gluon operators, "condensates", or (partial) conservation laws of 
(axial-) vector currents. These sum rules arc thus prime examples of connect- 
ing hadronic excitations to the underlying ground-state structure (symmetry- 
breaking pattern), including order parameters. We will briefly discuss two 
classes of sum rules. 
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2.3.1 Q CD Sum Rules 



QCD sum rules (QCDSRs) [55] are based on the analyticity of correlation 
functions, resulting in the dispersion relation 

r JmU^Uu) ^ ^^^ (6) 



The left-hand-side involves an integral over a (hadronic) spectral function in 
the time-like regime plus possible subtractions not indicated here. The right- 
hand-side is an expansion in space-like momenta, l/Q"^ {Q"^ = —q^ > 0), with 
leading perturbative terms and nonperturbative effects encoded in Wilson co- 
efficients, Cn, via quark and gluon condensates of dimension increasing with n. 
When applied to the light vector mesons, p and u, at finite temperature and/or 
density [56], it turns out that the largest sensitivity resides in the medium de- 
pendence of the four-quark condensates, {{qqY). Unfortunately, rather little 
is so far known about their temperature dependence from lattice QCD. The 
usual assumption is to factorize the four-quark condensates into a product 
of two-quark condensates corresponding to the assumption of "ground-state 
dominance" with an extra parameter, k, usually fixed in the vacuum, rep- 
resenting correlation effects. An application of QCDSRs to the p meson in 
cold nuclear matter is shown in Fig. 8 [58], indicating that the finite-density 
decrease of the condensates mandates a "softening" of the p spectral func- 
tion. The required low-mass enhancement can be satisfied by an increasing 
width, a decreasing mass, or a suitable combination thereof as indicated by 
the "allowed regions" enclosed by the bands in the right panel of Fig. 8, see 
Refs. [57,59,60]. 



2.3.2 Chiral Sum Rules 

Chiral sum rules (CSRs) [61,62] have been derived prior to QCD from current 
algebra and chiral Ward identities. The (partial) conservation of the (axial-) 
vector-isovector currents leads to relations between the pion decay constant, 
= 92 MeV (an order parameter of chiral symmetry breaking), to moments 
of differences between pertinent spectral functions. In vacuum they aT^^ 








^ The form of the sum rules in Eqs. (7)-(10) applies to the chiral Umit of vanishing 
current Hght-quark masses and pion mass. Corrections to the second Weinberg Sum 
Rule (WSR), Eq. (9), may not be smah [63,64,65]. The first WSR, Eq. (8), does not 
seem to be affected by finite quark mass corrections. 
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Fig. 8. QCD sum rule constraints on Breit-Wigner p spectral functions in the 
mass- width plane [58]. The bands indicate regions of mass and width values for 
which the deviation between the left-hand and right-hand sides of Eq. (6) is below 
0.2% (within the solid lines) or below 1% (within the dashed lines). The left column 
is for the vacuum case, which is used to determine the "correlation parameter" , n 
in the four-quark condensate while the right column is at nuclear-matter density 
(diamond: mass and width of the free p meson). 
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(8) 
(9) 
(10) 



where (r^) is the pion charge radius squared, Fa the axial vector form factor 
in the radiative pion decay, n Wij. The vacuum vector and axial vector 
spectral functions are 

= --Imn^5 . (11) 
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In Eq. (10), obtained in Ref. [66], {O) denotes a four-quark condensate. In 
the factorization approximation, {O) = (16/9) (gg)^. The direct connection of 
the CSRs to the vector correlator renders them particularly valuable in the 
context of dilepton production. The assessment of in-medium effects requires 
their extension to finite temperature, elaborated in Ref. [66]. Due to loss of 
Lorentz invariance when specifying the thermodynamic rest frame, the original 
vacuum results become energy sum rules at fixed three-momentum and split 
into longitudinal (L) and transverse (T) components of the correlators and 
quark condensate: 



/ 



Ap^(go,g)=0 , 



OO 

y'cigo'Ap^'^(go,g)=0, 



(12) 
(13) 



/ ql dql [Ap^(go, q) + 2 Ap^(5o, q)] = -47ra, + 2 {{0'% , (14) 



where Ap = pv — pA and (( ■ )) denotes an in-medium expectation value. The 
transverse and longitudinal components of the spectral functions are given in 
terms of standard projection operators. 



PV,A — Pv,A^T + Pv,A^L 1 



(15) 



where the pionic piece, p'^^ — f'^q^5{(f)Pi^ in vacuum, has been included in 
the longitudinal axial vector channel. In medium, the pion spectral function is 
subject to medium modifications as well. The in- medium Weinberg-type sum- 
rules (12)-(14) impose stringent constraints on both temperature and energy- 
momentum dependencies through the moments of chiral hadronic models for 
vector and axial vector spectral functions. The model-independent connection 
to lattice QCD can be implemented by employing the pertinent temperature 
dependencies of the pion decay constant and four-quark condensates which, 
in principle, are easier to compute than full spectral functions. 



We now turn to a more concrete discussion with examples of how to realize 
the three main points outlined at the beginning of this Section. 



17 




M„ [GeV] [GeV] M_,^ [QeV] 

Fig. 9. Comparison of electromagnetic emission spectra from QGP and HG at in- 



termediate Mil and qt in central Au+Au collisions at ^ys^ = 200 GeV. Left panel: 
IMR dileptons [67] including primordial Drell-Yan annihilation. Middle panel: direct 
photons [19] including primordial pQCD contributions. Right panel: low-mass dilep- 
tons with pair qt cuts. Note that the IMR dilepton spectrum (left) neither includes 
contributions from the (j) nor a cut on the sini^/e-electron rapidity. The single-elec- 
tron rapidity cut has been applied to the low mass region, causing a significant 
reduction in the yield. On the other hand, the low mass spectra (right) does not 
include dilepton emission from four-pion and higher states that dominate hadronic 
emission beyond the (p mass. 

2.4 Temperature and Degrees of Freedom 
2.4-1 Suitable Kinematic Regimes 

Thermal emission rates for electromagnetic radiation, Eqs. (1) and (2), can in 
principle be used to "infer" the temperature of thermalized matter in heavy- 
ion collisions if: (i) the emission strength represented by the EM correlator 
is reasonably well determined so that the T dependence essentially resides 
in the Bose factor and (ii) a kinematic window can be identified where ra- 
diation from a reasonably well-defined temperature regime prevails. Ideally, 
these conditions are met at the highest masses and energies where the cor- 
relators can be reliably evaluated in pQCD and thermal emission from the 
earliest phases dominates. In practice, however, the high-mass/energy region 
is dominated by dileptons and photons from primordial hard A^A^ collisions. 
At lower masses, the lower temperature contributions increase substantially. 
A more comprehensive discussion of the various sources will be given below. 

Fig. 9 shows an example of predictions for space-time integrated photon and 
dilepton spectra at intermediate masses and transverse momenta in central 
^/s^ = 200 GeV Au-|-Au collisions [19,67]. The thermal spectra are decom- 
posed into QGP and hadronic (HG) emission from an isentropically expand- 
ing thermal fireball assuming a critical temperature of = 180 MeV [68] 
where the total entropy is fixed to reproduce the observed hadron multiplic- 
ities at chemical freeze-out, /i^, Teh) = (25, 180) MeV. The pQCD photon 
rates, to leading order in as [69], and HTL-resummed dilepton rates [46] are 
convoluted over a chemically-equilibrated QGP assuming a formation time of 
To = 1/3 fm/c, translating into Tq = 370 MeV (670 MeV if the initial parton 
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densities are assumed to be undersaturated) . Uncertainties in the longitudinal 
expansion, affecting the QGP lifetime, can induce changes of the QGP spec- 
tra by up to 30%. The sensitivity to tq is larger, especially at masses above 
2 GeV [67]. The hadronic emission spectra include in-medium modifications of 
the EM correlator, see Sec. 2.5 for details, as well as chemical off-equilibrium 
in the hadronic evolution until thermal freeze-out, see Sec. 2.6. Three regimes 
emerge where QGP radiation outshines both HG and primordial emission, 
see Fig. 9): 

(i) three-momentum integrated dilepton spectra at intermediate mass, 1.5 < 
M < 3 GeV (left panel); 

(ii) direct photon spectra at intermediate transverse momentum, 1.5 < gj < 
3 GeV (middle panel); 

(iii) low-mass dilepton spectra at gt > 2 GeV ( "low-virtuality" photons; right 
panel). 

In practice a careful assessment of additional sources, including "pre-equilibrium" 
contributions and jet-plasma interactions such as Bremsstrahlung off quark 
jets or Compton scattering of gluon jets [70,71], is mandatory before firm con- 
clusions about the thermal component can be reached. We note that the pre- 
dictions for the thermal photon spectra, together with jet-plasma interactions 
and primordial photons extracted from pp collisions, compare favorably with 
preliminary RHIC direct photon data [72], as discussed later. Recent calcu- 
lations of the jet-plasma component suggest that it exceeds thermal emission 
for real and virtual photon transverse momenta g* > 3 — 4 GeV [70,71]. See 
Section 3.5 for a possible strategy for disentangling the different components. 

2.4-2 Direct Photons and Current RHIC Data 

To illustrate the uncertainties in and the required precision of a "tempera- 
ture measurement" , we describe an analysis of recent direct photon spectra in 
central 200 GeV Au-|-Au collisions at RHIC. The left panel of Fig. 10 com- 
piles several model calculations of thermal photon production mostly based 
on hydrodynamic evolution. The maximum initial temperature, Tq, is quoted: 
Srivastava et al. [73], tq ^ 0.2 fm/c and 450 < Tq < 660 MeV; Alam et 
al. [74], To = 0.5 fm/c, Tq = 300 MeVQ, Rasanen et al. [76], tq = 0.17 fm/c 
and To = 580 MeV; Turbide et al. [19], expanding fireball with tq = 0.33 fm/c 
and (To) = 370 MeV; d'Enterria and Peressounko [28], tq = 0.15 fm/c and 
To = 590 MeV; see also Steffen and Thoma [77], tq = 0.5 fm/c (To) = 300 
MeV. For similar initial conditions, the total thermal yields in these calcu- 
lations are compatible both with the data and each other within a factor of 

^ Alam et al. have recently [75] recomputed their hydrodynamic yields using higher 
initial temperatures, Tq = 400 MeV at tq = 0.2 fm/c, to improve the agreement 
with the data. 
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Fig. 10. Left panel: Thermal photon spectra for central Au+Au reactions at 
^/s^ = 200 GeV, computed within different models (see text), compared to the 
expected pQCD prompt 7 yields (T4^-scaled NLO pp calculations [9]: solid line, 
no symbols) and to the experimental total direct photon spectrum measured by 
PHENIX [8]. Right panel: Direct photon "nuclear modification factor", Eq. (16), 
for the 0-10% most central Au+Au reactions at ^/s^ = 200 GeV. The solid line is 
the ratio resulting from a hydro+pQCD model [28]. The points show the PHENIX 
data [8] over the same NLO yields while the dashed-dotted curves indicate the 
theoretical uncertainty on the NLO calculations. 

~ 2. While this confirms the dominant role of thermal radiation in the win- 
dow pt — 1-5 — 3 GeV, it shows that more quantitative, in-depth comparisons 
are required to disentangle the underlying assumptions on evolution model 
(boost-invariant hydrodynamics with or without transverse expansion, ther- 
mal fireballs, etc.) from production rates to narrow down the viable Tq range. 
The excess over that expected from primordial A^A^ collisions, the "prompt" 
contribution, is better illustrated by the nuclear modification factor, R\j^, the 
ratio of the direct photon spectra in AA collisions relative to either the pQCD 
prediction for pp collisions, scaled by the nuclear overlap interval, Taa, or, 
preferably, to pp spectra measured in the same experiment, 

07 /„ N dN^Iu/dpT . . 



shown on the right side of Fig. 10. The data are consistent with a significant 
excess over the next-to-leading order (NLO) pQCD predictions. However, we 
emphasize that below ~ 4 GeV it is not yet clear to what extent the NLO 
predictions in the denominator of Eq. (16) are applicabl^. In this regime, the 
calculated prompt yields are dominated by jet Bremsstrahlung, determined 
from the parametrized parton-to-photon GRV form factor [78], rather poorly 
known at the pertinent pt- The standard scale uncertainties in the NLO pQCD 

^ The denominator may be modified by the "isospin effect" or other factors, see 
Sec. 3.4. 
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calculations are ±20% for pr > 4 GeV but could become as large as ^200% 
for 1 < Pt < 4 GeV, as indicated by the dash-dotted lines on the right side of 
Fig. 10. Obviously, precise measurements of the direct photon baseline spectra 
in pp and d+Au collisions at ^/s = 200 GeV above pt = I GeV are essential 
for quantifying a thermal signal in Au+Au collisions. 

The measured slope of the resulting thermal photon spectrum, T^s, does not 
directly reflect the temperature of the hot matter since photons are emit- 
ted throughout the space-time volume of the evolving matter, implying vary- 
ing temperatures as well as blue shifts due to collective expansion. Never- 
theless, a correlation between the apparent photon slope and the maximum 
temperature attained in the system persists, as in the recent hydrodynamical 
study within [28]. The measured T^g provides an empirical link to the effec- 
tive number of degrees of freedom of the system via c?of = (30/7r^)(e/Tjj) or 
dof = (45/27r^)(s/Tjj) [28]. The initial maximum energy and entropy densi- 
ties are difficult to access experimentally. Indeed, all observables related to 
the initial energy and entropy densities such as the total transverse energy, 
the total particle multiplicity, and the colored-particle density encountered by 
quenched jets on their path through the medium are related to space-time 
averaged quantities. Information on the temperature dependence of doi can, 
however, be obtained via centrality and ^/s^ dependencies. It is possible to 
discriminate a QGP-like equation of state with a fixed dof above from a 
hadronic resonance gas with a rapidly rising number of degrees of freedom, 
by establishing the dependence of T^s on the charged particle pseudorapid- 
ity density [28]. As a further consistency check, one can relate do{ to suitable 
powers of the energy and entropy densities, dof oc s'^/e^ [27]. 

2.4-3 Chemical Ojf- Equilibrium 

While rapid thermalization of the matter at full RHIC energy is fairly well 
established, its composition in terms of quark and gluon degrees of freedom 
(chemical equilibration) is much less clear, as are the equilibration mecha- 
nisms themselve^. From the prevalence of gluons in the relevant x-range of 
the incoming nuclei at midrapidity, one might expect the early matter to be 
a gluon plasma (GP), as is routinely assumed e.g. in jet-quenching by radia- 
tive energy loss [79]. However, recent calculations of qq pair production within 
the classical fields generated by the incoming Au nuclei indicate a rather fast 
approach to chemical equilibrium [80]. This rapid equilibriation could have 
important consequences for disentangling the relevant fast parton energy loss 
mechanism (with a significantly reduced radiative loss for quarks due to the 
smaller color charge). EM probes are an obvious means of testing chemical 



' The expansion in hydrodynamic simulations is mostly driven by the ratio of pres- 
sure to energy density, P/e, where dof essentially drops out. 
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equilibration since gluons carry no electric charge [67,81,82,83,84]. LO pQCD 
processes (g + (7 — * 7 + X, qq ^ l^l^) in a thermally equilibrated but chem- 
ically off-equilibrium QGP suggest that the photon and dilepton production 
rates scale with A^^^ and A^ A^, respectively, where Aj are fugacities char- 
acterizing the deviation of the parton densities from chemical equilibrium 
(Ag = \q — \g = 1). At RHIC, a typical GP initial state with subsequent 
evolution using inelastic pQCD reaction rates starts from \g ~ 1/3, A^.^ < 0.1 
evolving to values of A ~ 0.5 or larger [84,85,86]. However, in an isentropic 
expansion with fixed initial entropy, undersaturated matter implies signifi- 
cantly higher initial temperatures at otherwise identical conditions; e.g. for 
central Au-|-Au collisions at RHIC, with tq = 1/3 fm/c, Tq ~ 370 MeV in a 
chemically-equilibrated system relative to Tq ~ 670 MeV in a system off chem- 
ical equilibrium. In the photon spectra, reduced fugacities in the emission rate 
are largely compensated by the higher temperatures in the QGP evolution 
with only a slight hardening of the slope parameter [26] . This effect appears 
to be more pronounced for thermal dilepton spectra, see the left-hand side of 
Fig. 9. Proving that a harder slope is evidence for a GP at RHIG is further 
complicated by the fact that the thermal yields are still fairly sensitive to the 
assumed thermalization time, Tq. A decrease in tq with chemical equilibration 
not only decreases the slope of the thermal spectrum but also increases the 
yield. Identifying a GP or more generally, the number of degrees of freedom 
above Tc with thermal dileptons will thus necessarily involve a quantitative 
assessment of both slope and absolute magnitude of the thermal spectrum, af- 
ter "removal" of non-thermal sources including Drell-Yan dileptons, correlated 
open-charm decays, as well as pre-equilibrium and jet-plasma interactions. In 
the following section we briefly summarize recent progress on the last two 
sources. 



2.4.4 Pre-Equilibrium and Jet-Plasma Emission 

Dilepton emission subsequent to the initial hard NN collisions but before the 
assumed Tq, the so-called pre-equilibrium contribution, can be addressed in 
a parton cascade approach. A corresponding calculation [89] predicted large 
emission rates, which, in fact, overestimate preliminary PHENIX data [90] 
(Fig. 11). However, if Landau-Pomeranchuk-Migdal (LPM) interference ef- 
fects are included, the yield is appreciably suppressed [87] , indicating that the 
combined initial pQGD plus pre-equilibrium yield is not very different from 
the pQCD contribution alone (Fig. 11), with the thermal yield dominating at 
Pt < 2.5 GeV, similar to the center of Fig. 9. Also note that the calculated 
thermal spectra in Fig. 11 agree reasonably well with those in the center panel 
of Fig. 9. However, the present model dependencies will have to be further re- 
duced to achieve enough sensitivity to discern the composition of the early 
matter and realize the desired temperature measurement. 
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Fig. 11. Comparison of (early) pho- 
ton emission from a parton cascade for 
central 200 GeV Au+Au collisions [87] 
with (squares) and without (diamond) 
LPM effects. Also shown are expanding 
fireball [88] and hydrodynamic [20] cal- 
culations of thermal photons, as well as 
preliminary data from PHENIX [90]. 



Fig. 12. Predictions for direct pho- 
ton spectra combining pQCD prompt 
photons, jet-plasma interactions [70] 
and thermal radiation [19] compared to 
prehminary PHENIX data [72]. (Here 
"Th-Th" and "HG" are thermal radia- 
tion from the QGP and hadron gas re- 
spectively.) 



Fig. 12 shows the combined prediction for direct photons from prompt (pQCD) 
production, jet-plasma interactions [70], and thermal QGP and HG radia- 
tion [19] (see center panel of Fig. 9) evaluated within the same expanding 
fireball. The comparison to preliminary low-p-r PHENIX data [72] is quite 
encouraging. In this calculation, jet-plasma contributions exceed the thermal 
yield at — 2 GeV, implying a rather narrow "QGP window". On the other 
hand, at — 5 GeV, where jet-plasma radiation [70] is still significant, the 
experimental direct-photon Raa in the right panel of Fig. 10 does not indi- 
cate much room for enhancement over prompt pQCD photons. Note, however, 
that thermal and jet-induced radiation are not independent contributions since 
both are affected by the lifetime and temperature or number and energy den- 
sity characterizing the QGP, imposing an additional consistency requirements. 
Furthermore, as emphasized in Refs. [70,71], the intensity of electromagnetic 
radiation from jet-plasma interactions is intimately related to radiative energy 
loss, i.e., jet quenching via gluon emission. In this context, the estimates for 
the jet-plasma photon and dilepton contribution in Refs. [70,71] constitute an 
upper limit, since high-p^ pion suppression has been entirely attributed to 
radiative energy loss whereas recent calculations find a significant role played 
by elastic energy loss [91], even within the same formalism [92]. 



2.5 In-Medium Spectral Functions below and above Tc 



In the mid 1990's, dilepton data from the CERN SPS triggered vigorous the- 
oretical activity in trying to assess modifications of vector-meson properties 
in hot/dense (hadronic) matter. The focus has been on the p meson due to 



23 



its prevalent role in dilepton emission, see, e.g. Refs. [29,30,31,32,33] for re- 
cent reviews. In Sees. 2.5.1 and 2.5.2, we briefly summarize some of the main 
features and insights that have emerged over the last ~ 10 years. 

2. 5. 1 Hadronic Many-Body Theory and Chiral Virial Expansion 

Effective hadronic models for vector mesons should be compatible with basic 
symmetry principles, most notably electromagnetic gauge invariance, vector- 
current conservation, and chiral symmetr^T^. In addition, it is essential that 
the underlying effective vertices are carefully constrained by phenomenological 
information such as hadronic and electromagnetic decay widths or scattering 
data. Measurements which provide information on in-medium effects near nu- 
clear matter density are particularly valuable, e.g. photoabsorption data on 
both nucleons and nuclei [93,94]. 

Most of the effective models with constraints built along the above lines have 
reached a reasonable degree of agreement with the data, generically predict- 
ing substantial broadening in matter with little mass shift0, characteristic 
of calculations both in cold nuclear matter [59,95,96,97,98] and in hot and 
dense matter [12,99], see Fig. 13 for two examples. Effective models also sug- 
gest that the effects of the baryonic component of the medium dominate over 
those from the mesonic one at comparable density, consistent with findings 
in IdJge-Nc QCD where meson-meson interactions are suppressed relative to 
meson-baryon ones). 

The broadening of the spectral functions, amounting to a total width of 
~500 MeV at nuclear matter saturation density and typically accompanied 
by a slight upward mass shift (left panel of Fig. 14), is in fair agreement with 
constraints from QCD sum rules [58], recall Fig. 8. When extrapolated to 
temperatures and densities close to the expected chiral transition, an almost 
complete "melting" of the p-resonance structure emerges. This is not only true 
for the net-baryon rich regime at SPS energies and below, but also in the cen- 
tral rapidity region at collider energies where the baryon chemical potential 
is small (center panel of Fig. 14). At the experimentally extracted chemical 
freezeout temperature, e.g. Tch — 180 MeV at RHIC, an appreciable density 
of baryon-antibaryon, BB., pairs is thermally excited [67] and mesons interact 
with both baryons and antibaryons. In addition, the notion of chemical freeze- 

^ In many instances, little is known about the chiral structure of baryonic and 
mesonic resonance couplings, especially if no pions are involved. 
^ A simple explanation of this feature is that imaginary parts of the in-medium 
self energies, which govern the broadening, are negative definite (ImS < 0) and 
therefore strictly sum up whereas real parts, which induce mass shifts, change sign 
around a resonance. Real parts therefore tend to cancel if the system is characterized 
by a rich excitation spectrum, as is the case for a hadronic resonance gas. 
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Fig. 13. Comparison of p spectral functions in hot hadronic matter within a many- 
body calculation [12] (gray lines, extending to M=0) and an approach based on 
imaginary parts of p-hadron scattering amplitudes supplemented by dispersion re- 
lations to obtain the real parts [99] (black lines). Reasonable agreement between 
the in-medium results is observed, especially for nucleon densities, n^v < 1 (in units 
of the nucleon saturation density, 0.16 fm~^). 

out implies that baryon-antibaryon annihilation in the subsequent hadronic 
evolution is suppressed [101]. Thus, the relevant quantity for medium effects 
on vector mesons is the sum of the B and B densities which, close to T^h, 
is quite comparable at = and Hb = 250 MeVH^. The baryon-density 
effects on the p are most pronounced at masses below ~ 0.5 GeV instead of at 
and above the free p mass (compare the long-dashed and short-dashed lines 
in the center panel of Fig. 14). While the (j) appears to be less sensitive to 
the baryonic component of the medium, this conclusion may be altered once a 
better understanding of recent photon- and proton-nucleus (f) production data 
has been achieved [102,103]. 

The p and possible u "melting" has interesting implications that deserve fur- 
ther theoretical investigations: 

(i) The very short mean-free path of the p (and other hadrons) close to is 
suggestive of hadronic liquid formation [104]. Thus from a hadron structure 
objective, the matter properties in the phase transition region may change 
rather smoothly from a hadronic liquid to the sQGP liquid. 

(ii) The hadronic in-medium EM correlator, based on a "melted" p, is surpris- 




Experimentally, the total baryon rapidity density, dN^_^§{y = 0)/dy, is indeed 
very comparable at the maximum SPS {y/s^ = 17.3 GeV) and RHIC i^/s^ = 
200 GeV) energies. The total hadron rapidity density (mostly due to pions) is a 
factor of ~ 2 larger at RHIC, implying an accordingly lower total baryon density 
at the transition. However, most of the pertinent medium effects on the p spectral 
function build up at densities at or below qq, see the left panel of Fig. 14. 
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Fig. 14. In-medium vector meson spectral functions within hadronic many-body 
theory. Left: the p in cold nuclear matter [100] at several densities proportional to 
the saturation density, qq = 0.16 fm~^. Center: p behavior at RHIC conditions close 
to Tc with and without baryon effects [67]. Right: same as the center panel but for 
uj and (j) [67]. 

ingly similar in shape and magnitude to the QGP correlator evaluated in 
HTL pQCD [46] at all masses [12,105]. This is also suggestive of a rather 
continuous transition from HG to QGP close to Tc, even on the level of 
spectral functions. The approximate coincidence of the bottom-up and top- 
down extrapolated hadronic and pQCD calculations, while not understood 
theoretically, enhances our confidence in applying these emission rates to 
fireball and hydrodynamic evolution across T^, necessary for the space-time 
integrated thermal dilepton spectra, rendering the dilepton spectra rather 
insensitive to the exact value of Tc [106]. 

Similar conclusions about the p also emerge from the chiral virial approach [94] 
where medium effects on the vacuum vector correlator are evaluated within 
a pion- and nucleon-density expansion coupled with vacuum Vtt and VN 
scattering amplitudes constrained by chiral symmetry. The p peak is quenched, 
though not broadened, and its low-mass shoulder is substantially enhanced, 
predominantly due to baryon effects in the heat bath [107]. This agreement, 
at least at low and moderate densities and temperatures, is a consequence of 
the constraints imposed on the underlying hadronic interactions. 

2.5.2 Dropping Mass 

Models involving dropping vector-meson masses [108] have recently been revis- 
ited within the so-called vector manifestation of chiral symmetry [32,109,110]. 
Using the Hidden Local Symmetry (HLS) framework, where the p mass is 
generated via a Higgs mechanism, an alternative representation of the chiral 
group has been proposed in which the chiral partner of the pion is the lon- 
gitudinal p meson, the so-called "vector manifestation" of chiral symmetry 
instead of the conventional realization where the a is the chiral partner of 
the pion. The HLS approach results in a satisfactory vacuum phenomenology. 
A renormalization group analysis with hadronic loop effects reveals a fixed 
point with a vanishing vector coupling constant. When applied to the second 
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order finite temperature chiral phase transition, matching the vector and axial 
vector correlators to the operator product expansion (space-like q^) requires 
that bare p mass vanish at the critical point, becoming degenerate with the 
pion mass, which persists when carried on-shell due to the fixed-point nature 
of the transition. The vector dominance model (VDM), which works well in 
the vacuum, is perhaps violated at finite temperature [110], suppressing tttt 
annihilation to dileptons via an intermediate p, replaced by direct annihila- 
tion via intermediate photons. The violation of vector dominance could make 
the observation of a dropping p mass difficult in dilepton spectra, perhaps 
reconciling the new NA60 data [6,7] with a dropping-mass scenario. Ref. [Ill] 
calculates the dilepton rates in the vector manifestation scenario. Despite vi- 
olation of the VDM, the rates clearly exhibit a p peak with dropping mass, 
at least up to T ~ 0.85Tc. We emphasize that the matching procedure is 
only vahd sufficiently close to Tc and that "ordinary" hadronic medium effects 
become dominant at lower temperatures. In the low-temperature limit, mod- 
ification of the p mass to leading order in temperature, 0(T^), is at variance 
with chiral symmetry. Thus, the notion of a "flash temperature" has been 
introduced [33,112], below which the "intrinsic" temperature dependencies of 
the parameters in the Lagrangian are void. 

At the quark- ant iquark level, interactions in the vector channel are believed 
to be rather weak since the p mass is nearly twice the constituent quark 
mass. Finite temperature effects resulting in enhanced interactions and an 
accordingly reduced p mass are not easily conceived [113]. (See Ref. [114] 
for an alternative view.) It would also be interesting to work out how the 
presence of hadronic many-body effects (especially baryons), as discussed in 
the previous section, affect the matching procedure and the resulting (axial) 
vector spectral functions. 



2.5.3 Resonances in the sQGP 

Another interesting development perhaps related to EM measurements are 
the conjectured hadronic bound states in the (s)QGP [35,36]. Detecting sig- 
natures of vector states above Tc in the dilepton spectrum hinges on whether 
their mass is sufficiently large, MyiT > T^) > 1 GeV. As elaborated after 
Eq. (3), QGP radiation can only compete with or dominate contributions from 
the longer-lived and larger-volume hadronic phase at large masses, especially if 
the resonance structure depends on temperature. The IQCD spectral functions 
and dilepton rates shown in Figs. 5 and 7, respectively, indeed indicate reso- 
nances with masses in the M ~ 2 GeV ~ lOTc regime, roughly scaling with 
temperature. The existence of these states may be understood [35] to be due 
to heavy-quark quasiparticles bound by a rather strong color Coulomb-type 
attraction, also present in heavy quarkonium states. In this case, heavy-quark 
symmetry implies approximate degeneracy of vector ("p" or J/ip) and pseu- 
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Fig. 15. Ratio of dilepton spectra at 
RHIC from a resonance-enhanced QGP 
rate to perturbative qq annihilation [37] 
obtained from the QGP and mixed 
phases of an expanding thermal fire- 
ball [67,68]. The first peak is the vec- 
tor mass close to Tc while the second 
is in the vicinity of the "zero-binding" 
line at T ~ 1.5 — 2Tc. The x-axis is in 
units of the quark-quasiparticle mass, 
Mq ~ 1 GeV. The quark-quasiparticle 
width is assumed to be 0.1 GeV [36]. 



Dilepton Excitation Function in Central Au-Au (N ^^.^=330) 




Fig. 16. Excitation function of low-mass 
dilepton spectra for central Au-|-Au col- 
lisions from top SPS to top RHIC en- 
ergy. The upper (lower) set of four 
curves corresponds to thermal radia- 
tion from hadronic matter (QGP) us- 
ing in-medium p, oo and (j) spectral func- 
tions [12,68] (HTL annihilation [46]), 
assuming a critical temperature be- 
tween 175 MeV at ^/s^ =17 GeV and 
180 MeV at > 100 GeV. 



doscalar ( "vr" or rjc) states. The connection between these resonances and the 
pion mass dropping to near zero when approaching Tc from above has been 
addressed in Ref. [114]. 

Quantitative signatures of vector resonances above in the dilepton spectrum 
have been evaluated in Ref. [37] . Convoluting the temperature-dependent res- 
onance decays with an expanding fireball at RHIC [67,68] predicts an enhance- 
ment over baseline pQCD emission scenario {qq annihilation) by about a factor 
of ~ 2, see Fig. 15. This result is sensitive to the vector-meson width which in 
turn is governed by the width of the quark quasiparticles. Smaller widths lead 
to narrower peaks and thus a stronger enhancement over the pQCD spec- 
trum. The quasiparticle width is expected to be ~ 0.2 GeV based on self- 
consistent solutions of a qq scattering equation [36] using input interaction 
potentials from finite temperature IQCD or Nambu-Jona-Lasinio four-quark 
interactions [115]. The corresponding elastic scattering rates of ~ l/(fm/c) 
are suggestive of the short thermalization times deduced from hydrodynamic 
analyses of elliptic flow measurements and could therefore provide a link to 
the early thermalization puzzle at RHIC. 
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2.6 Low-Mass Dilepton Spectra 

To illustrate the predictions of medium- modified vector mesons within effective 
models, and in particular to investigate the importance of the baryonic compo- 
nent of the medium, we summarize a pertinent excitation function of low-mass 
dilepton spectra in central Au+Au collisions in Fig. 16. Thermal dilepton rates 
in the QGP [46] and HG [12,67] phases are convoluted with isentropic fireball 
evolution similar to that underlying Figs. 9 and 15, assuming a chemically 
equilibrated QGP which converts into a chemically equilibrated HG at {iJ>%, 
Tc) values compatible with: (i) thermal models for hadron production in cen- 
tral AA collisions [25] and (ii) a total entropy that reproduces available par- 
ticle multiplicity data [2,3,4]. The hadronic evolution subsequent to chemical 
freeze-out is augmented by effective chemical potentials for hadrons that are 
stable under strong interactions {e.g. tt, K, rj, baryons and antibaryons) . This 
is mandatory to maintain the observed chemical composition until thermal 
freeze-out [116] and, in particular, implies sizable total {B + B) baryon den- 
sities in the later stages of the hadronic evolution (recall Sec. 2.5.1). Surpris- 
ingly, there is no large change in either shape or magnitude of the (hadronic) 
dilepton-spectrum excitation function for 20 < -y/s^ < 200 GeV. The main 
reasons for this outcome are: 

(i) Despite the large range in baryochemical potentials (yU^ = 25 — 250 McV 
at chemical freezeout with = 175 — 180 MeV) and thus in net baryon 
density, the prevalent baryon-induced medium effects are comparable once 
the B + B density is properly accounted for. 

(ii) The lifetime of the hadronic (and mixed) phase changes little since the 
larger volume expansion at higher energies is essentially compensated by an 
increase in radial flow inherited from the QGP phase. 

This scenario should be contrasted with one where the medium effects are 
sensitive to the net baryon density such as simple dropping mass parameter- 
izations, my/mv = (1 — Cqb/qo) x (1 — {T/Tc)'^Y^^. In this case, a stronger 
variation of the excitation function is anticipated, with weaker effects at higher 
collision energy. QGP emission increases appreciably with ^/s but remains sub- 
dominant (<20%) in the low-mass region at all energies if no significant qt cut 
is applied, as may be expected for larger initial temperatures and longer QGP 
lifetimes. 

At masses above ~1 GeV, the hadronic EM spectral function is dominated 
by four pion and higher contributions, encompassing annihilation reactions 
such as pp, nuj, and nai. While these contributions are not included in the 
hadronic matter calculation shown in Fig. 16, they may become significant 
at M > 0.9 GeV where a similar enhancement could be related to effects of 
partial chiral symmetry restoration [117,118]. We will return to this issue in 
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Sec. 2.7.2. 



Finally, a few remarks on models of the space-time evolution are in order 
since these provide crucial input on the thermodynamic parameters for the 
equilibrium EM emission rates. Hydrodynamic models, if applicable, are the 
approach of choice since they are formulated using the same variables as the 
thermal emission rates. The present RHIC data suggest that ideal hydrody- 
namics gives a good approximation of the first ~ 5 fm/c after thermalization, 
encompassing the QGP, "mixed" phases and possibly the early hot-|-dense 
hadronic liquid phase for T > 150 MeV. For lower temperatures, however, vis- 
cosity effects are expected to become significant [119]. It is presently not clear 
how finite viscosity affects calculations of dilepton and photon emission. The 
underlying uncertainties must be scrutinized, especially since the low-mass and 
low-momentum spectra receive significant contributions from later stages. One 
can either implement viscosity effects in hydrodynamics to retain the notion 
of thermodynamic variables or switch to transport theory [120,121,122]. In the 
latter option, it is nontrivial to properly implement broad resonances [123]. Al- 
ternatively, local temperatures and (baryon) densities could be extracted from 
transport simulations and convoluted with the equilibrium EM emission rates. 
The relative agreement of this method with viscous hydrodynamics could pro- 
vide an estimate of the uncertainties in the integrated EM spectra and used 
to better calibrate fireball models, which are suitable parametrizations of mi- 
croscopic evolution. Such comparisons will become particularly relevant if less 
penetrating probes, e.g. tttt or 777 invariant-mass spectra, are calculated, as 
discussed later. Transport-based approaches are mandatory for thermal freeze- 
out to quantitatively account for finite sizes, lifetimes and mean-free paths. 
In the late 1990s, the agreement of hydrodynamic [124]. transport [125] and 
fireball models [12,126] with CERES low-mass dileptons [10] at the SPS has 
been reasonable, albeit with somewhat limited theoretical and experimental 
precision. 



2. 7 Chiral Symmetry Restoration 
2. 7. 1 Direct and Indirect Approach 

We now address the question of how, in principle, in-medium effects detected 
in dilepton spectra can be used to draw conclusions about chiral symmetry 
restoration (xSR). An unambiguous consequence of the xSR is that isovector, 
vector and axial vector correlation functions, which are very different in the 
vacuum, become degenerate at and beyond the chiral transition. The ques- 
tion is how this happens, see Fig. 17 for an illustration. We reiterate that 
the effects of chiral symmetry breaking are concentrated at low masses since, 
already in vacuum, the correlators become degenerate in the pQCD regime. 



30 



0.08 



0.06 



g 0.04 

> 

C 0.02 




-> 2mi V 1 

p(770) +cont. 
aj(1260) + cont. 







s [GeV ] 



3 

o 

Oh 

00 



a 

3 

o 



Dropping Masses ? 



pert. QCD 



Mass 



Melting Resonances ? 




Mass 



Fig. 17. Left: isovector, vector and axial vector spectral functions in the vacuum as 
measured in hadronic r decays [127] with model fits using p and oi resonances plus 
a perturbative continuum [128]. Right: two schematic scenarios of chiral symmetry 
restoration in hot and dense matter. 

Therefore, symmetry breaking constitutes an inherently nonperturbative phe- 
nomenon which, ideally, can be addressed with input from experiment, theo- 
retical models and lattice QCD computations. 



A direct way to search for xSR is measurement of the in-medium axial vec- 
tor spectral function in connection with model comparisons, as also done for 
dileptons. Measurements of the 7r^7 invariant- mass spectra were suggested 
[129] as a means of searching for xSR- This was partly motivated by similar 
measurements of tt+tt" invariant- mass spectra in heavy-ion collisions which 
indicated p modifications in the late stages of peripheral 200 GeV Au+Au 
collisions at RHIC [130]. Absorption effects on the outgoing pions limit the 
TTTT information to rather dilute stages while ttj spectra may probe somewhat 
deeper into the fireball. However, emission from later collision stages is ad- 
vantageous since a narrower density and temperature window is probed so 
that the convolution over the space-time history becomes less of an issue. Ex- 
perimentally the challenges are the rather low rateJ"! and the rather broad 
structure of the ai(1260) resonance (~ 0.4 GeV in vacuum), making it suscep- 
tible to distortions due to background subtractions. Simulations will be shown 
in Sec. 4.2.2. 

In addition to the direct experimental approach to the axial vector chan- 
nel, we now formulate a well-defined theoretical procedure based on chiral 
hadronic models to connect experimental information on the vector correla- 

Even though the radiative decay branching ratio of the oi is small, the (not so 
small) absolute decay width, Tai^n-y ~ 0.7 MeV, is the relevant quantity for thermal 
radiation. 
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tor (dileptons) to (first-principle) information on xSR from lattice QCD. The 
in-medium versions of the chiral sum rules, Eqs. (12), (13) and (14), critical 
to the calculation, are obtained as follows: 

(1) First, calculate the vector (V) and axial vector (A) spectral functions as a 
function of temperature and density in a chirally invariant model, including 
as many of the constraints as possible (see Sec. 2.5.1). 

(2) Then insert the spectral functions into the Weinberg sum rules to evaluate 
the temperature dependence of the pion decay constant and four-quark con- 
densate and compare to lattice QCD results. Note that f-K{T) and {{qqY){T) 
are presumably more easily evaluated in IQCD than spectral functions since 
IQCD is primarily applicable to the finite temperature axis, i.e., at /Ug=0, 
the closest relation between IQCD and heavy-ion experiments is realized in 
the central rapidity regions at RHIC and the LHC. 

(3) Finally, perform detailed comparisons of the in-medium effects on the vector 
correlator with dilepton data as a function of centrality, ^/s^, mass and qt- 
spectra. The comparison requires additional input from realistic expansion 
models {e.g. hydrodynamical and transport simulations), which have been 
thoroughly tested against the large body of hadronic observables. 

The three different energy moments of py — probed by the chiral sum 
rules, provide detailed constraints on the energy dependence of the in-medium 
spectral functions. In addition, each in-medium chiral sum is valid for a given 
three-momentum, providing further kinematic information. Therefore, if a chi- 
ral hadronic model complies with both theoretical (2) and experimental (3) 
tests, a tight connection between lattice QCD and data has been established, 
producing explicit evidence for chiral symmetry restoration without a direct 
measurement of the axial vector correlator. In the absence of (unquenched) 
lattice data for (low-mass) spectral functions for at least the next ten years, 
a systematic approach involving effective models is the only way to interpret 
data in terms of x^R- Experimental guidance is crucial for progress in under- 
standing the underlying nonperturbative physics. 



2.7.2 LMR-IMR Transition: Chiral Mixing 

In a low-temperature pion gas, the expectation values of vector and axial vec- 
tor correlators can be evaluated modcl-indcpcndcntly based on chiral reduc- 
tion formulae in connection with a low-density expansion. The leading medium 
effect has first been derived in Ref. [131] in the chiral limit (m^ = m^^d = 0) 
and amounts to chiral correlator mixing. 



nv(g) = (l-e)n™^(g) + enr(g) 




(17) 
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Fig. 18. Comparison of NA60 dimuon data [6,7] in central In (158 AGeV) + In 
collisions with calculations. Left: an expanding fireball calculation [117,118] for QGP 
and hadronic emission. The hadronic emission is based on in-medium p [12], uj [67] 
and (j) spectral functions as well as four-pion contributions with an enhancement 
from chiral mixing according to Eq. (17). Right: hydrodymanic convolution [132] of 
the dilepton rates following from the chiral virial expansion [94,107]. 

where e = encodes the thermal pion density (smaller for > 0). 

Interactions of the vector current with pions from the heat bath quench the 
vacuum vector correlaton I and its admixture of the axial vector correlator 
induced hj V + n ^ A and A -\- n ^ V processes. The axial vector current 
is analogously affected. Full mixing corresponds to e = 1/2, implying degen- 
erate correlators. Broadening and possible mass shifts of the p are thus due 
to higher-order effects in both T and g^. The mixing has the interesting fea- 
ture of fil ling in the "dip" in the s = 1 — 2 GeV^ region of the vacuum vector 
correlatoillfj (see the left panel of Fig. 17). Dilepton enhancement by up to 
a factor of 2 over the vacuum vector spectral function in this mass region is 
therefore a signature of the approach to chiral restoration, via ttoi annihi- 
lation (four-pion contributions not present in the vacuum EM correlator) or 
QGP emission. While the earlier SPS dilepton data did not have the neces- 
sary precision for the required quantitative analysis (recall Fig. 2), it may be 
feasible with the new NA60 data [6,7] shown in the left panel of Fig. 18. The 
blue dashed curve, a theoretical upper estimate [118], employs Eq. (17) with 
e(T) = in7r(T)/?T,^(Tc) where n^(T) is the pion density at T < Tc, including 
pion chemical potentials below = 175 MeV, and removing the ai —>■ irp 
decay, included in the p in-medium spectral function. On one hand, it is grati- 
fying to see that this calculation properly accounts for the excess spectrum in 
the relevant mixing regime. On the other hand, the data are still reasonably 
well described without mixing. To be sensitive to the mixing effect, both the 
data and the theory require an accuracy of at least 20%. The chiral virial ap- 



The shape of the correlator is unaffected. 

Full mixing in this regime leads to degenerate V and A correlators that closely 
coincide with the pQCD qq continuum level, interpreted as lowering the "duality 
scale" from s ~ 2.5 GeV^ in the vacuum to about 1 GeV^ for full mixing. It is further 
tempting to interpret the p melting found in hadronic many-body theory (recall 
Sec. 2.5.1) as lowering the duality scale as s ^ 0, implying chiral restoration [105]. 
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Fig. 19. Effective slope parameters as a function of dimuon invariant mass for excess 
radiation in In (158 AGeV) + In collisions at the SPS. Note the different vertical 
and horizontal scales on each plot. Left: experimental results with dNch/dy > 30 
summed over all centralities and fit over 0.4 < < 1. 8 GeV/c [135]. Center: 
fireball calculations of semicentral collisions based on thermal emission plus free 
p decays at freezeout using the local slope at qt = 1.5 GeV [136]. Right: fireball 
calculations of semicentral collisions based on thermal and non-thermal sources fit 
over 1.0 < qt < 1.8 GeV. The lower dashed line is a baseline calculations while 
the upper dashed and dot-dashed lines also include t-channel meson exchange in 
irp TT/i/i reactions [106]. 

proach, when folded over hydrodynamic evolution, as shown in the right panel 
of Fig. 18, also describes the region M > 1 GeV well [132]. The free EM cor- 
relator with mixing effects is a key ingredient, resulting in an enhancement 
consistent with maximal mixing as the upper limit and the free EM correlator 
as the lower limit. Although the p peak is quenched (Sec. 2.5.1), the abscence 
of p resonance broadening results in a ~ 40% overestimate of the yield around 
the free p mass. The enhancement below the p mass is again accounted for, 
with important baryon contributions. In both approaches underlying Fig. 18, 
the QGP yield is small. The conclusions of Ref. [118] on p broadening and 
the importance of baryon-driven medium effects in the context of the NA60 
data have been confirmed in Ref. [133] which, however, attributes the bulk of 
the enhancement above the p mass to QGP radiation. The relative QGP to 
four-pion yields in the IMR is essentially determined by the choice of Tc in 
the fireball evolution [106]: QGP dominates the IMR if = 160 MeV while 
hadronic contributions are dominant if Tc = 175 MeV. Irrespective of whether 
the source is of QGP or hadronic origin, the IMR enhancement is associated 
with matter at temperatures close to Tc. The robustness of this conclusion is 
again a consequence of "parton-hadron duality" in the underlying emission 
rates. 

Additional information on the nature of the emitting source might be obtained 
from quantitative analysis of dilepton qt spectra, including elliptic flow [134]. 
Such analysis has become possible with the NA60 data [135]. Effective inverse 
slope parameters, Tg^, have been extracted from dilepton excess spectra in the 
momentum range 0.4 < < 1.8 GeV, obtained from inclusive In-|-In collisions 
at the SPS with dNch/dy > 30. The resulting Tes are shown as a function of 
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dilepton invariant mass on the left side of Fig. 19. The extracted values of 
Teff appear quite large relative to the hadronic slopes at freezeout, requiring a 
surprisingly strong collective flow for the moderate system size at SPS energies. 
The maximum Teff(M), in the free p mass region, is indicative of p decays in 
the late stages of the collision where the line shape is expected to approach the 
vacuum shape while the effective temperature, T^s — Tfo + m-u^, is primarily 
due to the blue shift imprinted by the collective expansion velocity, v. Note 
that the slope of the excess radiation at M = rrip ,^ is larger than that of the uj 
(lower circle at M = mp ,j, suggesting that the p freezes out significantly later 
than the r/, oj and (the upper open circle at M = m^^^, estimated for p decays 
around freezeout). In the IMR, Tefj is significantly reduced, in agreement with 
thermal emission early in the evolution, from temperatures close to Tc where 
the flow has not yet developed much, implying a small blue shift. Quantitative 
theoretical analyses [106,136] are rather involved, see the center and right 
panels of Fig. 19, especially for momenta above qt — 1 — 1.5 GeV where 
nonthermal sources are expected to become significant [137], e.g., Drell-Yan 
dileptons or primordial p decays. 

It would be very valuable to obtain similar information for Au+Au collisions at 
RHIC. The thermal yield is expected dominate further, contributing at higher 



momenta, due to the larger system size and collision energy I. Unfortunately, 
the increased contribution from correlated charm decays is sensitive not only 
to the total charm cross section but also to the charm momentum spectra and 
thus to charm thermalization (at low and intermediate pt) and energy loss 
(at high Pt)- However, open heavy flavor spectra are interesting independent 
probes of the medium and their measurement therefore constitutes one of the 
central goals of RHIC-II, see the heavy flavor part of this report. It would also 
be very illuminating to check whether the much increased partonic collectivity 
at RHIC reflects itself by an increased in the IMR. 
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2.8 Electromagnetic Signatures of the Color Glass Condensate 



Saturation physics has been applied to the description of RHIC data quite 
successfully, from hadron multiplicities and the phenomenon of limiting frag- 
mentation in Au+Au and d+Au collisions to the produced hadron transverse- 
momentum spectra in d+Au collisions at mid- and forward rapidity (for recent 
reviews and an extensive list of references, see [138]). Nevertheless, in order 
to establish gluon saturation as the dominant physics responsible for these 
phenomena at RHIC and beyond, and to rule out other phenomenological 
scenarios, one needs to consider further tests of the CGC formalism such as 

The multiplicity, dNch/dy in Au+Au collisions at RHIC is about a factor of four 
larger than that of In+In collisions at the SPS. 

We thank Nu Xu for an interesting discussion on this point. 
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Fig. 20. Dilepton production in quark-nucleus scattering. 

the predictions of saturation physics for electromagnetic processes. In this sec- 
tion we outhne electromagnetic signatures of the CGC at RHIC. Specifically, 
photon and dilepton production in d+Au collisions are considered. These pro- 
cesses have an advantage over hadronic processes in that they do not interact 
strongly after they are produced and the nonperturbative hadronization pro- 
cess is absent. Furthermore, photon and dilepton production can shed light 
on the validity of the recombination model approach to hadron production 
which are also capable of fitting the available data, albeit with a few as- 
sumptions. Since saturation physics predicts similar suppression patterns for 
photon, dilepton and hadron production in d+Au collisions, an experimental 
confirmation of this generic prediction would be a major step in establishing 
saturation physics at RHIC and in ruling out recombination as the physics of 
hadron production in d+Au collisions. 



2. 8. 1 Dilepton and Photon Production 

We consider the dilepton-production cross section in quark-nucleus scatter- 
ing [139] 



q{p) + q{q) + + {k^) + X , (18) 

shown in Fig. 20, where ki and /c2 are the momenta of the two leptons. 
Photon production can be obtained by taking the photon virtuality (dilepton 
invariant-mass) to zero. The differential cross section at fixed impact param- 
eter, ht, is given by 



d%t d^ktdlnM'^ dz 



37r 



d^q 



zHl 



-zil - z) 



+ M2(l - Z)][{kt -Zltf + M2(l - Z)] 
1 1 



n 2 



[kj + M2(l - z)] [{h -zltY + M2(l - z)] 



(19) 



where z is the fraction of the incoming quark light cone energy carried away 
by the (virtual) photon while xa is the Bjorken x probed in the target nucleus. 
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All quark masses arc ignored, is the dilepton invariant mass squared with 
h — (lt + kt and kt is the transverse momentum of the lepton pair. 

Eq. (19) is the standard expression obtained if propagation of the incident 
quark through the nucleus is assumed to be eikonal, i.e., the transverse mo- 
mentum transferred to the nucleus by the incident quark is much less than its 
longitudinal momentum [140]. All the information on the degrees of freedom 
in the target is contained in the dipole cross section, a^^^^i^{bt,lt,XA). Gluon 
saturation physics comes in via the dipole cross section, which is determined 
by the evolution equations of the Color Glass Condensate. Given an initial 
condition for the dipole cross section, the JIMWLK equations in the Color 
Glass Condensate formalism determine the dependence of the dipole cross 
section on the collision energy (or alternatively, xa) and It- This is the main 
difference between the CGC formahsm and the other approaches [140]: while 
the CGC formalism can predict the xa and It dependence of the dipole cross 
section, other approaches cannot and must instead motivate a suitable form 
from phenomenological considerations. We note that to obtain the invariant 
cross section for photon production, da'^'^^'^'^'^ /(fbtcfhtdz, from Eq. (19), 
the mass is set to zero, M — 0, and the dilepton vertex factor, a/STrM^, is 
removed. 

Since the dipole cross section is also the main ingredient in the hadron pro- 
duction cross sections in deuteron (proton)-nucleus collisions, a similar nuclear 
modification factor, RdA, is predicted for dilepton production and hadron pro- 
duction. Parametrically, one expects the nuclear modification factor, RpA, to 
scale with the nucleon number A like RpA ~ where 1 — 70 ~ 0.628 is the 

BFKL anomalous dimension. It should be emphasized that in QCD, small- 
X evolution (BFKL) is the only way to generate leading-twist shadowing if 
partonic degrees of freedom (quarks and gluons) are used. 

Dilepton production has an additional knob to turn in order to change the 
kinematics and probe QCD dynamics in different settings, the dilepton in- 
variant mass. The nuclear modification factor, calculated using a saturation 
inspired model of the dipole cross section [141], is shown in Fig. 21. 

Since the dilepton production rates are small due to the electromagnetic cou- 
pling, we consider the invariant-mass dependence of the transverse-momentum 
integrated cross section. The integration over Eq. (19) can be done analytically 
to obtain 



(j^(jqA-*ql+l- X 

' d^dM^dz 



1 — z 



37r2 ^2 



[1 + (1 - zf] KH^^^Mtt] + 2(1 - z) Kl[^^^^ 



Mrq 



(20) 
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Fig. 21. Nuclear modification factor, RdA, of dilepton production in quark-nucleus 
scattering [141]. The dotted, short-dashed and long-dashed curves are for M = 2 
GeV and y = 0.5, 1.5 and 3 respectively. The solid curve is for M = 4 GeV and 
y = 3. 

To relate this to proton (deuteron)-nucleus scattering, Eq. (20) must be con- 
voluted with the quark (and antiquark) distributions, g(x, M^) (g(x, M^)) in 
a proton or deuteron. This can be written in terms of the proton (deuteron) 
structure function F2 [140], 



daP 



A^l+l- X 



a 



1 



(PhdM'^dxF Qtt'^ Xq + x„ 



dz I drl Fl{xJ z) (jl,^^^^{xg, b, tt) 
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(22) 



and 



F|^^e5x[g;(x,M2) + g/a;,M2)] 
/ 



is the proton structure function. 

The invariant dilepton production cross section is shown in Fig. 22 for proton- 
proton, proton-nucleus and deuteron-nucleus collisions as a function of dilep- 
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Fig. 22. Dilepton production cross sec- 
tion at y = 2.2. 



Fig. 23. The dilepton nuclear modifica- 
tion factor in pA and dA collisions at 
y = 2.2. 



ton invariant mass, while Fig. 23 shows the nuclear modification factor for 
both pA and dA collisions at RHIC. The calculations are for fixed rapidity, 
y = 2.2, and the most central collisions [142]. The HKM parameterization of 
the deuteron wave function is used for the deuteron projectile, a few percent 
effect. 

Photon production and photon and hadron correlation functions have been 
calculated using the same formalism [142]. The hadron-photon cross section 
is a very sensitive probe of the dipole profile, the main ingredient of single- 
particle production cross section in the CGC formalism. Therefore, experimen- 
tal studies of electromagnetic probes in deuteron-nucleus collisions at RHIC 
can shed light both on the dynamics of gluon saturation and on the role of 
saturation physics in the observed suppression of the hadron spectra in the 
forward rapidity region at RHIC. This measurement can also clarify the role 
of hadron recombination models at RHIC, at least in the forward rapidity 
region, since recombination effects will not be present in electromagnetic fi- 
nal states. Therefore, observation of dilepton or photon suppression would be 
strong evidence for the CGC. 

It is worth noting that, since the saturation scale of the proton (deuteron) is 
very small at mid- and forward rapidity at RHIC, CGC predictions for proton- 
proton collisions will have large uncertainties. Indeed, as recently shown [143], 
particle production in proton-proton collisions cannot be reliably calculated at 
RHIC, since the saturation momentum is small, while d+Au calculations are 
under much better quantitative control. Therefore, in any saturation-inspired 
calculation of the nuclear modification factor, the pp cross section in the de- 
nominator should be understood to be a fit to the data while the d+Au cross 
section can be calculated using saturation physics and used to determine the 
physics dominant in the forward-rapidity region at RHIC. This, in turn, will 
have significant ramifications for the LHC. 
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Fig. 25. Dilepton mass spectrum from 
PHENIX [146] compared to a cocktail of 
hadronic sources and charm decays. 



3 OBSERVABLES: STATE OF THE ART 



3.1 Low-Mass Dileptons 



The p meson, with its short hfetime (1.3 fm/c in vacuum), large dilepton 
decay width (Tee = 7 keV) and prevalent coupling to tttt (and possibly qq) an- 
nihilation has long been identified as the most promising probe of in-medium 
modifications of hadron properties close to the QGP phase boundary. Changes 
in p mass or width were anticipated as precursors of the chiral transition more 
than two decades ago [144]. Dielectron mass spectra from CERES [10,22] over 
the last ten years have shown a significant (factor of 2 — 3) excess over known 
hadronic sources in the mass region below the free p mass - a result which 
spurred vigorous theoretical activity. However, the CERES experimental un- 
certainties, including limited statistics and mass resolution, could neither con- 



firm nor refute scenarios with medium effects I. Recently, the NA60 experi- 
ment at the SPS, a significant upgrade of NA50 including high-precision track- 
ing close to the vertex in a high magnetic field, measured dimuon mass spectra 
in 158 AGeV In+In collisions with an unprecedented ~ 20 MeV mass resolu- 
tion and high statistics in both the low and intermediate mass regions [6,7]. 
The results in the light vector-meson region are shown in Fig. 24. After sub- 
tracting the "hadron decay cocktail" with t], uj, and decays, the excess 



A significant improvement in mass resolution was achieved with the CERES Time 
Projection Chamber upgrade [145]. 
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spectrum has a clear albeit rather broad peak around the nominal p mass, 
as shown in Fig. 18 and compared to theoretical predictions using in-medium 
p spectral functions with either a dropping mass [11] or a strong broaden- 
ing, as calculated in hadronic many-body theory [12]. The excellent resolution 
and statistics, made possible by detector upgrades and high luminosity, can 
now clearly distinguish between these two approaches and be considered as 
a benchmark for future dilepton measurements. (How well other approaches 
can explain the measured data remains to be seen.) Measurements of simi- 
lar precision - made possible by RHIC-II - will be essential to further improve 
our understanding and elevate it from the qualitative to the quantitative level, 
necessary for drawing conclusions about the nature of chiral symmetry restora- 
tion. 

Comparable quality measurements of the p spectral function do not yet exist 
at RHIC. PHENIX is measuring dielectrons in the central arm. The mass 
spectrum in minimum bias 200 GeV Au-|-Au collisions is shown in Fig. 25 
along with the hadronic cocktail and the background from charm. While the 
J/ip measurement is very clean, the resolution in the low mass region is not 
adequate to test in-medium modifications of vector mesons. Such precision is, 
however, expected once the hadron-blind detector, one of the major PHENIX 
upgrades becomes operational (see Sec. 4.1.2). 

STAR measured the 7r"^7r~ invariant- mass spectra [130] and found indications 
that the p mass is both px and multiplicity-dependent in pp, d+Au and periph- 
eral Au+Au collisions. At low transverse momenta, 0.5 < px < 1-5 GeV/c, the 
p mass peak is 3 — 8% lower. Whether this is caused by a p mass shift [147,148] 
or other effects {e.g. Bose-Einstein correlations, thermal phase space in con- 
nection with resonance broadening, underlying "cr" decays, etc. [129,149]) still 
needs to be clarified. The measurement was insensitive to possible changes in 
the p width [150] which requires an improved determination of the background 
sources including the "physics background"! . 

The combination of leptonic and hadronic p decays is a valuable tool for 
disentangling in-medium modifications due to the hot and dense phases from 
the more dilute phases. The leptonic decays probe the entire collision history 
while the hadronic decays probe the fireball surface and the late stages. 

3.2 The (f)(1020): Hadronic vs Leptonic Decays and 

Since the mass is barely twice the kaon mass, — 2mK — 30 MeV, a 
study of hadronic, KK, relative to leptonic, e+e~, decays has long 

In a different energy range, KEK (E325) measured a drop of the p mass wliereas 
JLAB (CLAS) reported broadening without any drop of the p mass. 
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Fig. 27. Elliptic flow as a function 
of PT for and A in 200 GeV 

Au+Au collisions [153]. 



been suggested as a sensitive test of a dropping mass or braodening of the 
(j) width, possibly related to partial chiral symmetry restoration. If the mass 
drops, the dominant (p — ^ KK channel is quickly kinematically suppressed and 
the ratio of production cross sections of the hadronic, B{(j) — > KK) = 49.2%, 
and leptonic, 5(0 e^e~) = 3 x 10~^, channels, (t(0 — > KK)/a{(f) e'^e"), 
should change dramatically. Similar arguments may apply if a reduced kaon 
mass or inelastic scatterings such a.s (p + n —>■ K + K* lead to a strong in- 
crease of the width. In this case, regeneration has to be accounted for as 
well [151]. Recent preliminary results from PHENIX [152] suggest that, when 
normalized by the branching ratio, dN{(f) e~^e~)/dy is larger and the rriT 
slopes are steeper in AA relative to pp collisions. However, the experimen- 
tal uncertainties prevent definite conclusions. The main problem is the small 
signal-to-background ratio (S/B) in the dielectron channel which can only be 
cured by better background rejection. However, the yields from —>■ KK 
extracted by STAR [153] are systematically higher than those measured by 
PHENIX in the same decay channel. Careful analysis shows that pt spectra 
measured by the two experiments are not inconsistent [154]. 

A difference was observed between —>■ KK measured by NA49 and — >■ ee 
measured by NA50. Recent CERES [155] results for both decay channels in 
Pb+Au collisions and measurements of ^ f^~^f^~ decays in In+In collisions 
from NA60 [156] are consistent with data from NA49. Thus the difference was 
not due to the observed decay channel but the method of observation. This 
experience emphasizes both the importance of powerful detector upgrades and 
the ability to measure different decay channels in the same apparatus over a 
wide momentum range. 
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Fig. 28. Preliminary PHENIX direct 
photon spectra in 200 GeV Au+Au 
collisions, for all centralities and min- 
imum bias collisions. The curves are 
NLO pQCD calculations with the scale 
set equal to pT- 
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Fig. 29. Preliminary PHENIX direct 
photon spectra in 62 GeV Au+Au colli- 
sions at different centralities compared 
to NLO pQCD predictions at scales of 
0.5pT, Pt, and 2.0pT- 



Since the (p meson mass is comparable to some baryon masses, it can also 
play a major role in distinguishing between the mass and species dependence 
of quantities such as the nuclear modification factor, Rcp, and elliptic flow. 
Recent results from STAR [153] in 200 GeV Au-|-Au collisions show that both 
Rcp (Fig- 26) and V2 (Fig. 27) are consistent with parton recombination and 
collective flow at the partonic level, as well as with the PHENIX data [157]. 



3.3 Direct Photon Spectra 



The latest direct photon spectra from Au-fAu collisions over the full centrality 
range are compared to NLO pQCD calculations at 200 GeV (Fig. 28) and 
62GeV (Fig. 29). While the data seem to favor a perturbative scale of 0.5pr, 
particularly at lower transverse momenta, they are consistent with a scale 



equal to px at higher momenta I. One explanation could be the onset of 
additional photon sources such as photons from jet-medium interactions at 
intermediate px (see Fig. 12). Before drawing any conclusions, high quality 
pp data have to determine the level of agreement of the data with the NLO 
pQCD predictions, see Sec 3.4. Unfortunately there are different competing 
sources of the excess relative to baseline pQCD calculations proposed - all with 
their respective theoretical uncertainties. Studying inclusive cross sections is 
insufficient for disentangling these sources. A possible solution, involving much 
larger data samples than currently available, is proposed in Sec. 3.5. 
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A similar trend is observed in the n spectra [158]. 
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Fig. 30. Direct photons in 200 GeV 
pp collisions from different PHENIX 
runs divided by a fit to the preliminary 
Run- 5 results. The black circles show 
preliminary Run-5 results to illustrate 
the quality of the fit while the red cir- 
cles are published Run-3 data. The blue 
curves are NLO pQCD calculations. 
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Fig. 31. Prehminary PHENIX di- 
rect photon Raa in the most central 
(0-10%) 200 GeV Au+Au collisions, fit 
to the pp direct photon spectra mea- 
sured in the same experiment. 



The first RHIC measurement of the photon nuclear modification factor, shown 
in Fig. 1, suggests that the direct photon Raa is neither enhanced nor sup- 
pressed, as opposed to the observed strong hadron suppression shown on the 
right side of Fig. 1. The vr" suppression observed in the very first RHIC Au+Au 
run but absent in d+Au collisions was interpreted as a final-state effect. Both 
the TT*^ and direct photon Raa interpretations assume that, at high where 
hard scattering should dominate, scaling the pp cross sections by the nuclear 



overlap function T^_bl£J is a sensible estimate of the expected AA yields in the 
absence of a medium. Jets are quenched relative to this expectation. Strictly 
speaking, only the direct photon Raa proves that this expectation is reason- 
able: photons produced by hard scattering should leave the medium unchanged 
and their yield should thus scale with Tab- This assumption has indeed been 
confirmed within ~ 20 — 30% uncertainty by RHIC Run-2 Au+Au data [8] 
and by the integrated Raa with pt > Q GeV/c where Raa was calculated using 
an NLO pQCD calculation since pp direct photon data were not yet available. 

The much larger Au+Au data set from Run-4, the two different collision en- 
ergies (a/s/vtv = 200 and 62 GeV) and measurement of pp reference data in 
the same experiment with similar systematic errors made it possible to study 
the evolution of Raa with p^ and collision centrality in much greater detail. 
The fit to the pp data differs both in shape and magnitude from NLO pQCD 
[159]. The data are higher, particularly at the lower and upper ends of the 
5 < Pt < 23 GeV/c range, see Fig. 30. The Raa obtained from a fit to 



Colloquially, but quite misleadingly often called "scaling with the number of bi- 
nary collisions" or A'cou scaling 
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Fig. 32. Preliminary PHENIX direct photon Raa in the most central (0-10%) 62 
GeV Au+Au collisions, using an NLO pQCD calculation with scales of 0.5pT, Pt, 
and 2.0pT- 

the pp data is shown in Fig. 31 for the most central 200 GeV Au+Au colli- 
sions. There is an apparent suppression at very high px- At least part of the 
suppression has been predicted due to the "isospin effect" [160], a natural con- 
sequence of the different quark content {a oc Se^) of protons and neutrons "^^ . 
Hence the suppression in the photon Raa is referred to as "apparent" . Other 
mechanisms such as modifications of the parton densities and fragmentation 
functions which enhance or deplete high pt direct photon yields in heavy-ion 
collisions have also been proposed. Finally, distinguishing direct photons from 
tt" decay photons producing an overlapping shower in the PHENIX calorimeter 
becomes increasingly difficult at Pt > 14 GeV/c. 

However, RHlC's ability to provide collisions over a wide energy range may 
allow some means of separation. If the suppression at 200 GeV, as shown in 
Fig. 31, is indeed only an isospin effect, it should scale with xt = Iprj^- 
Therefore, in lower energy collisions it should manifest itself at lower p^ where 
the experimental difficulties due to overlapping showers are absent. PHENIX 
measured a similar photon Raa in 62 GeV Au+Au collisions using NLO pQCD 
calculations for pp (preliminary data from PHENIX on direct photon produc- 
tion in pp at 62GeV will be available soon.) The results for the 10% most 
central collisions is shown in Fig. 32. While within errors it is consistent with 
unity, using measured data in the ratio may change this conclusion, as was 
the case at 200 GeV. 

Although not an electromagnetic probe per se, the vr" Raa in Figs. 33 and 34 

'^^ This implies that while pp data are the proper reference for the hadron Raa, in the 
photon case the proper mixture of Upp, (7p„, and (T„„ must be used. For example, in 
minimum bias collisions <7^a/-^co11 = {^/A'^){Z'^app + 2Z{A — Z)apn + (A — Z)'^ann)- 
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Fig. 33. Preliminary PHENIX tt^ Raa 
in central 62 GeV Au+Au collisions us- 
ing a fit to the world average of previous 
pp data (black circles) and a PHENIX 
pp measurement (red circles). 
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Fig. 34. Preliminary PHENIX 7r° Raa 
in central 62 GeV (red circles) and 200 
GeV (orange triangles) Au-|-Au colli- 
sions. Both ratios use PHENIX mea- 
surements. 



0.25 


' 1 ' 


0-20% 








20-40% 




0.2 
0.15 






























0.1 
















0.05 




























^— 















% H 






-0.05 




















40-60% 






0-92% 




0.25 
0.2 












b.g.) 




0.15 












inclusive 7 


) 


0.1 

0.05 














" 


{ IT 












-0.05 
















1 






o R v,(i 


nclusive y) 


- Vj(b.g.) 





2 4 60 2 4 6 



PT(GeV/c) 

Fig. 35. The measured V2 of inclusive photons (^2^'^ solid circles) and expected 
photon V2 from hadronic decays (^2^, open squares) [163]. The subtracted quan- 
tity, ^ — v^^ is plotted at the bottom of each panel (open circles). Here 
R = (A^'dir 7 + Abg)/Abg corresponds to the product fg"^ {R — 1)). Four centralities 
are shown, including the minimum bias result. 



illustrates the importance of measuring the reference in the same experiment 
with similar systematic errors. Fig. 33 compares the previous and the current 
preliminary 7r° suppression in 62 GeV Au+Au collisions. The earlier result 
used a fit to the ISR pp data (upper, black circles) whereas the new result uses 
the pp cross sections now measured in the same experiment (PHENIX). The 
two results in Fig. 33 would clearly lead to very different physics conclusions. 
Fig. 34 shows the new 62 GeV result compared to the Raa obtained in 200 
GeV collisions. 



46 



3.5 Photon Azimuthal Asymmetries (Elliptic Flow) 

Fig. 1 showed that the overall azimuthally-integrated high-p^ photon yields 
scale with A^cou ^-nd are well described by pQCD within current experimen- 
tal and theoretical uncertainties. However, this global agreement may mask 
more subtle effects. It is even possible that the agreement is only accidental, 
due to cancellations of processes that enhance and quench the photon yield. 
An important step toward clarifying the level of agreement is the study of 
azimuthal asymmetries in the photon distributions, specifically their elliptic 
flow, V2- If photons from the initial hard scattering do not interact with the 
medium, their V2 is expected to be zero. However, initial hard scattering is 
not the only source of photons in AA collisions. They may also originate from 
jet-thermal parton interactions or from Brehmsstrahlung off a quark. These 
photons are expected to exhibit a negative V2 [161,162] (since more material 
is traversed out-of-plane, the major axis in coordinate space) with a strong 
PT-dependence. However, photons from thermal radiation should reflect the 
dynamical evolution of the hot and dense matter, thus carrying a positive V2. 

The first measurement of photon elliptic flow at RHIC is shown in Fig. 35 
[163]. The measurement is quite delicate due to the large background from 
tt" decay-photons that inherit the V2 of the parent 7r°. While the measured 
V2 of inclusive photons is consistent with the V2 of hadronic decay photons, 
i.e. a zero net direct photon flow, the errors are appreciable and the direct- 
to- inclusive photon ratio is very small at low pr- The quality of the current 
data is insufficient to verify the predictions [161,162], including the sign of 
the net flow. Much higher statistics can help remedy the situation, at least 
at higher px- Although the net direct photon flow is predicted to decrease, 
the statistical errors will become smaller and, equally important, the direct- 
to-inclusive photon ratio will increase dramatically. But even at high px the 
net flow will be a competition between processes with positive and negative 
V2. New analysis techniques may be able to statistically disentangle isolated 
and non-isolated direct photons in heavy-ion collisions. Jet-photon conversions 
primarily produce isolated photons [161] with negative V2. The magnitude of 
this flow depends strongly both on px and the jet energy-loss mechanism in 
heavy-ion collisions. Therefore, a measurement of the isolated photon V2 may 
provide an independent constraint on energy-loss models. 



3.6 Electron Raa and Flow 

The recent measurement of "non-photonic" single-electron spectra, associated 
with decays of open heavy-flavor hadrons, at RHIC lead to two unexpected 
and very important results: (i) the nuclear modification factor, Raa, shows a 
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Fig. 36. Left: illustration of the direct photon measurement via low- mass dilep- 
tons [72]. Right: preliminary results on low- and intermediate-pT direct photons, 
the "photon excess ratio" , with the traditional calorimeter method (red circles) and 
low mass electron pairs (blue circles). 



strong suppression in central Au+Au collisions [164,165,166], comparable to 
the suppression observed for pions; (ii) the elliptic flow, V2, is signiflcant at 
low pt, up to 10%. Although such strong suppresion was predicted a decade 
ago [167], it was later argued that the suppression would be mild due to the 
"dead cone effect" for heavy quarks. Measurements from both STAR and 
PHENIX have, in fact, shown strong charm suppression and signiflcant charm 
flow (early thermalization) in 200 GeV Au+Au collisions: observations not 
simultaneously explained by current theories. These issues primarily concern 
charm and bottom physics and are discussed in great detail in the heavy flavor 
part of this Report. However, they may be signiflcant for "classic" electromag- 
netic probes as well since heavy flavor suppression actually aids measurements 
of intermediate mass continuum dilepton radiation by reducing the combina- 
torial background from charrj^ in this regime, see Fig. 39. 
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3. 7 Direct Thermal Photons via Low-Mass Dielectrons 



A very promising approach for measuring low-p^ direct photons is to utihze 
low-mass electron pairs from "internal conversions", as was first applied in 
heavy- ion collisions by PHENIX [72]. The basic idea is that any process pro- 
ducing a real photon can also produce a very low mass virtual photon [169] 
which subsequently decays into an e^e~ pair. This direct photon signal com- 
petes, of course, with dielectrons from vr", 77, etc Dalitz decays. The dielectron 
rates and mass distributions are described both for low-mass direct photons 
and Dalitz decays by the Kroll-Wada formula [169], 



Ami f 2ml\ 1 , , 9,,,/ ml\^ , , 

1-^1 + ^ ^^ee l-irf , 23 



1 C^A^ee _ 2« 

drriee 37r ^ 



where the form factor, F, is unity for real photons. Note that the phase space 
for Dalitz decays is limited by the mass of the parent meson, m^.^ < M^rO 
while for direct photons it is not, ~ Pt- Therefore, the measurement 
becomes relatively clean for pj- > 1 GeV/c, still in the low-pj' realm where 
"traditional" calorimeter measurements have serious difficulties. The method 
is illustrated in the left panel of Fig. 36 while the resulting photon excess ratios 
are shown on the right. The systematic errors are much smaller for the new 



method I. The statistical errors become large for > 4 GeV/ c although ~ 15 
times more data have been analyzed for dielectrons (solid red circles) than for 
the calorimeter analysis (open blue circles). There is an important lesson here: 
even signals such as thermal photon radiation, typically not considered to be 
"starving for statistics", novel but promising new analyses techniques require 
large luminosity increases. 

Before drawing any conclusions about the thermal nature of the apparent 
direct photon signal at low pt the baseline in pp has to be established. This is 
shown on the left panel of Fig. 37 (final results from PHENIX): the signal is in 
good agreement with NLO pQCD predictions. In contrast, in Au+Au (right 
panel) there is clearly an additional source at low px- 

The conversion approach is also appealing from the theoretical point of view 
since the framework outlined in Sect. 2.1 accommodates a comprehensive 
treatment of real and virtual photons, as represented by Eqs. (1) and (2). 
However, additional theoretical work is needed for a more complete assess- 
ment of possible very- low mass dielectron sources [21]. 



This background is irreducible for the PHENIX HBD. 

The errors are smaller, at least for the direct photon excess ratio. However, the 
absolute normalization and the direct photon cross sections are difficult to assess 
and therefore have been inferred using the calorimeter. 
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Fig. 37. The fraction of the direct photon component as a function of px in 200GeV 
pp (left) and Au+Au (right) colHsions. The error bars and the error band are the sta- 
tistical and systematic uncertainties, respectively. The curves are from NLO pQCD 
calculations. 



4 THE FUTURE 



4-1 Detector upgrades 



The original designs of both PHENIX and STAR anticipated later upgrades to 
enhance their physics capabilities. In fact "upgrading" as a response to spec- 
tacular rises in luminosity and physics insights from data already taken was 
and remains an almost continuous process from the very start of operations. 
Few can be tied uniquely to RHIC-II with the exception of upgrades facili- 
tating the handling of high luminosities and data rates. Nonetheless, several 
major projects should be discussed here because of their magnitude and their 
impact on the future capabilities of RHIC-II, see Table 1. These projects are 
in different stages of development. Some are in the early R&D phase while 
others are ready to be installed soon. Here we concentrate on upgrades that 
are particularly important for electromagnetic probes. The STAR Time-of- 
Flight (TOF) detector and Heavy Flavor Tracker (HFT) improve the elec- 
tron identification and the displaced vertex measurement by rejection of elec- 
trons from charm, bottom and Dalitz decays as well as conversion electrons. 
The PHENIX Hadron Blind Detector (HBD) will reject Dalitz pairs for light 
vector meson measurements. The Silicon Vertex Detector (SVTX) will mea- 
sure displaced vertices to identify electrons from open charm. The Nose-Cone 
Calorimeter (NCC) will measure photons and vr^'s at forward rapidities. A 
new high-resolution sampling calorimeter is proposed for STAR to measure 



50 



direct photon correlations. In the following, we will discuss these systems in 
more detail. 



4.1.1 STAR Time- of -Flight and Heavy Flavor Tracker 

STAR electron identification is made possible by a combination of energy loss, 
dE/dx, by charged particles due to ionization of the time-projection chamber 
(TPC) gas and a velocity measurement with the TOP system [170]. The rela- 
tivistic rise of the electron dE/dx separates electrons from hadrons except at 
the crossovers with pions at ~ 0.2 GeV/c, kaons at ~ 0.6 GeV/c, protons at 
~ 1.1 GeV/c and deuterons at ~ 1.5 GeV/c. A time-of-flight measurement, 
with the requirement that |1 — /3| < 0.03, eliminates slow hadrons and cleans 
up the crossovers, resulting in clean electron identification [170]. 

In addition to direct measurements of open-charm hadrons via Ktt decays, 
the STAR HFT [171] will be a powerful means of discriminating primordial 
electrons from background electrons in the measurement of electromagnetic 
probes. 

As discussed earlier, direct photon and lepton production is rare and over- 
whelmed by photons and leptons from electromagnetic decays of hadrons and 
7 conversions to electrons, 7 — > e~^e~. The 7 conversions, to a large extent, 
occur in the detector material. The HFT detector will reduce the background 
electrons and positrons from these 7 conversions. By requiring hits in the HFT, 
electrons from photon conversion outside the HFT, i.e., in the Silicon Strip 
Detector (the upgraded STAR Silicon Vertex Tracker) and the TPC inner field 
cage, are rejected. To estimate the signal-to-background ratio in the vector- 
meson measurements, STAR adopted a reasonably conservative approach and 
assume that the HFT can reject 7 conversions by a factor of 10 [171]. Al- 
though the configuration of HFT has evolved significantly, the simulation of 
background rejection is not expected to change. Another electron background 
arises from semileptonic decays of heavy quarks. If heavy-quark spectra are ex- 
trapolated from pp collisions, the dominant dilepton source in central Au-|-Au 
collisions at intermediate mass [68] is due to semileptonic cc decays. With the 
large charm yield at RHIC, the latter are comparable to the yield from 7 con- 
versions and 7r° and f] Dalitz decays after HFT rejection. Detailed simulations 
show that the HFT is capable of rejecting ~ 75% of e~^e~ pairs from decay 
while preserving 50% of the direct e+e~ pairs. 

The large reduction in electron background will enable STAR to observe 

the electromagnetic signal from low-mass vector mesons and radiation of 
intermediate-mass dileptons with a few hundred thousand central Au-|-Au 
events. The rejection of 7r° and rj Dalitz decays by a factor of three (single 
track) can be achieved by measuring both electrons of the pair, made possible 
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Fig. 38. Dielectron invariant mass distributions for central 200 GeV Au+Au colli- 
sions. The solid black curve is the prediction for thermal emission [68] in the STAR 
acceptance. The red curve at the top is the total dielectron invariant-mass spectrum 
in the 2004 configuration assuming electron PID from full TOF coverage. The pink 
dashed line is the e^e~ pair spectrum from semileptonic charm decays, the grey 
curve is the charm e'^e~ spectrum after HFT DCA < 80 /um, and the dot-dashed 
line are from vr'^ and r/ Dalitz decays after rejection by the TPC. 

by the large acceptance of the STAR TPC. Future work on further rejection of 
such dilepton-background pairs with the inner tracker alone will be exploited. 
With the upgrades, STAR expects to detect 6K and 22K u decays in 200 
million recorded central Au+Au collisions. These are to be compared with 
~ lOK for the and ~ 6K for the uj in central In-|-In collisions, presented by 
NA60 [156]. 

Fig. 38 summarizes the dielectron background and signal invariant-mass dis- 
tributions. The signals of medium-modified vector mesons and thermal QGP 
radiation (black curve) are from calculations [68] folded over the STAR ac- 
ceptance. The uppermost (red) curve is the total dielectron invariant-mass 
spectrum in the 2004 configuration, obtained from the single-inclusive elec- 
tron spectrum measured in 200 GeV Au+Au collisions assuming electron PID 
from full TOF coverage. The pink dashed line is the e^e~ spectrum from 
semileptonic charm decays derived from the non-photonic single-electron spec- 
tra measured in Au+Au collisions. The gray curve is the charm e^e~ distribu- 
tion after applying the HFT distance of closest approach cut, DCA < 80 /im. 
The dot-dashed line is 7r° and r] Dalitz decays after TPC rejection. The net 
result is a signal-to-background ratio that, even in the continuum regime, is 
< 0.1, comparable to the central In+In NA60 measurements [7]. The standard 
way of dealing with the residual background is the mixed-event method, used 
by CERES and NA50/NA60 at the SPS and will be used in PHENIX and 
STAR. 

With the proposed Time-of-Flight, Heavy-Flavor Tracker and Data Acquisi- 
tion System (DAQIOOO) upgrades, STAR will be able to take data at a rate 
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e*e' invariant mass specta, Au+Au, 200 GeV, dN|,„/dy = 650 




Fig. 39. Combinatorial background for low-mass electron pairs compared to the to- 
tal signal from vector mesons and charm with and without the HBD. Blue squares: 
total e^e~ combinatorial (HIJING), no rejection from the HBD. Red circles: combi- 
natorial background after rejection with HBD. Blue stars: combinatorial background 
from charm alone (PYTHIA), not rejected by the HBD, measured separately by the 
Silicon Vertex Detector (SVTX). Solid red line: total combinatorial background af- 
ter HBD rejection (HIJING-I- PYTHIA). Dashed blue line: pure signal from direct 
vector-meson and Dalitz decays after freezeout "cocktail" without the charm contri- 
bution (EXODUS). Solid black line: total signal from vector mesons, Dalitz decays 
and charm (EXODUS-^PYTHIA). 

of 1000 Hz with very little dead time. At the same time, the collision vertex 
has to be limited to ~ ±5 cm due to the acceptance of the HFT [171]. The 
current data-taking rate is 50 Hz with a collision diamond of ~ ±50 cm. With- 
out a machine upgrade, the average luminosity delivered to STAR is 8 x 10^^ 
cm~^s~^ or 6 kHz Au+Au minbias nucleus collision rate. Taking into account 
the factor of 3 — 5 loss from the vertex constraint and centrality binning 
{e.g. 10% most central), central-triggered Au+Au events can be recorded at a 
rate of about 200 Hz with all the available luminosity, significantly below the 
DAQIOOO capability. A factor of 10 luminosity upgrade will enhance both the 
collider and detector. In peripheral collisions, the statistics will be lower due 
to lower event multiplicities. However, triggering will be more effective even 
for low-pT lepton pairs. Since the luminosity decreases quadratically with de- 
creasing beam energy, an efficient energy scan at RHIC requires a luminosity 
upgrade, as discussed in more detail in Sect. 4.3. 



4.1.2 PHENIX Hadron Blind Detector and Silicon Vertex Detector 

As discussed earlier, the measurement of the low-mass electron pair contin- 
uum to infer the in-medium modifications of the light vector-meson (p, u, and 
(f)) spectral functions is crucial for understanding the fate of hadron masses 
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and the approach to chiral symmetry restoration in the hot, dense matter 
created in heavy-ion colhsions. However, the e+e~ signal is overwhelmed by 
the combinatorics from 7r° Dalitz decays and photon conversions. The combi- 
natorial background from open-charm decays is also significant, see Fig. 39. 
Fortunately, the typical opening angle of the e^e~ pair is very small, both for 
7r° Dalitz decays and photon conversions. Thus unless they are swept apart 
by a magnetic field, a proximity cut on electron pairs is an effective veto on 
Dalitz and conversion electrons. The PHENIX magnet has an inner and an 
outer coil which are typically powered by currents in the same direction for 
maximum field to optimize tracking to the highest possible pt- However, the 
current in the inner coil can be reversed, making the magnetic field vanish- 
ingly small around the collision region, allowing the installation of a Hadron 
Blind Detector (HBD) to eliminate the combinatorial background from Dalitz 
decays and conversions. This detector, described in Ref. [172], is a windowless 
Cherenkov detector operated with pure CF4 in proximity focus configuration 
with a Csl photocathode and a triple GEM detector with pad readout. It 
is a crucial upgrade for the exploration of the properties of the new matter 
created. As illustrated in Fig. 39 it will reduce the combinatorial background 
and increase the signal/background ratio from ~ 10~^ to ~ 10~^. The latter is 
again very comparable to NA60 and the expected performance of STAR after 
upgrades, described in the previous section. 

The remaining background from open charm will be separately measured in 
PHENIX by a Sihcon Vertex Detector (SVTX), which will measure the heavy- 
fiavor displaced vertex with a 40 //m resolution of the distance of closest ap- 
proach. The SVTX resolution is driven by the cr of 123 //m and 462 ixm for 
and decays, respectively. The SVTX will have a central barrel and two 
endcap detectors, covering both central and forward rapidities and providing 
inner tracking with full azimuthal coverage and up to \ri\ < 2.4. This, in partic- 
ular, will enable the measurement of correlated e/i invariant-mass spectra and 
thus provide a standalone determination of the correlated open heavy fiavor 
component of the dilepton spectra. 



4.1.3 The PHENIX Nose- Cone Calorimeter 

PHENIX can measure photons, vr^'s and ?7's well at midrapidity but the rapid- 
ity coverage of the electromagnetic calorimeter in the central arm is limited 
to \y\ < 0.35, making full jet reconstruction very difficult. Also, several large 
y measurements in d-l-Au collisions suggest that \ow-x gluons might be satu- 
rated so that the CGC describes the results, including hadron suppression at 
large rapidity relative to no suppression at y = 0. If the large y suppression 
is indeed a consequence of initial-state gluon saturation, photons should also 
be suppressed, see Sec. 2.8. Such a study is only feasible with a calorimeter at 
large rapidities. The limited acceptance of the central arm also makes crucial 



54 



7— jet measurements very difficult. 

PHENIX tlius proposed to replace the current copper magnet nosecones by 
a calorimeter covering 1 < I77I < 3. This Nose-Cone Calorimeter (NCC) is 
a silicon-tungsten sandwich sampling calorimeter, longitudinally segmented 
in three sections (two electromagnetic and one hadronic) and read out by 
1.5 X 1.5 cm^ Si pads. In addition, two layers of 468 fim pitch strip-pixels 
are added to achieve separation of direct photons from high-energy 7r°'s by 
shower-shape reconstruction. The total depth is ~ 35Xq radiation and 1.3Labs 
nuclear absorption lengths. The expected resolution for electromagnetic show- 
ers is 23%/ V^. The longitudinal segmentation allows to distinguish between 
electromagnetic and hadronic showers. The coverage is sufficient to reconstruct 
the entire jet energy. 

Therefore, jet physics and energy-loss studies using both photon-tagged jets 
and leading tt'^'s will be possible far away from central rapidity with the NCC. 
It will also be very useful for studying P-state quarkonia via Xc ~^ jJ/i' 
and possibly Xb 7^ radiative decays together with the muon spectrometer, 
giving information on the different quarkonium dissociation temperatures and 
feed down contributions to the J/ip and T yields [173,174]. The NCC is also an 
important part of the spin program by measuring the polarized gluon structure 
function at low x. 



4.2 New Measurements 

4.2.1 Direct Photon Correlations (HBT) 

The HBT correlation of photons from 7r° or rj decays is exceedingly small 
because of the "large" distances, > 10'' fm, at which the decay occurs. Typ- 
ical correlations from 7r° decay photons are /xeV/c in relative photon four- 
momentum. Therefore, any measurable HBT correlation will come from pho- 
tons directly emitted from the collision. 

Photon HBT correlations encode space-time information about the system 
emitting the photons. Many calculations, as well as "simple" quantum me- 
chanics, indicate that by studying HBT in different pt regions, emitting sys- 
tems of different temperatures and thus at different stages of evolution can be 
identified [175]. In addition, interferometry can also be used to measure the 
yield of direct photons [18]. 

While the effect of the 7r° decay background is to dilute the HBT signal, the 
decay in itself does not generate a fake correlation. However, the residual cor- 
relation of decay photons due to HBT 7r° correlations could cause a fake HBT 
correlation. When the measured HBT parameters determined from charged 



55 



pions are added into a Monte Carlo calculation of the residual correlation, 
the change is negligible compared to reasonable estimates of the direct photon 
correlation assuming that the charged and neutral pion HBT parameters are 
the same [176]. 

The dilution effect, on the other hand, is a serious issue. Reasonable estimates 
of the visible HBT correlation indicate that A, the effective amplitude of the 

HBT correlation, is only a few parts per thousand. Consequently, great care 
must be taken in designing experiments to measure the direct photon HBT. 
Corresponding work on this problem is briefly summarized in the remainder 
of this section, based on a more complete, technical paper [176]. 

STAR studies showed that the best way to measure the HBT correlation is 
to use one photon that converted in a thin (10% Xq) converter placed inside 
the inner-field cage of the TPC and another (real) photon detected in the 
barrel calorimeter. It is thus possible to be sensitive to photons with nearly 
identical directions. However, in order to make useful measurements of the 
HBT parameters, the energy resolution of the electromagnetic calorimeter 
must be substantially improved over the resolution of the present calorime- 
ter, SE/E = 15% /VE. Pertinent studies show that 6E/E = 3%/VE is 
certainly sufficient and some preliminary Monte Carlo studies indicate that 
even 5%/y/E may be adequate albeit near the limit of usefulness. Work is 
in progress to determine the precise requirements, along with studies of the 
"Shashlyk" design, a relatively cheap calorimeter where alternating plates of 
radiator and scintillator are read out by an array of wavelength-shifting fibers 
which run parallel to the long axis of the calorimeter and penetrate the ra- 
diators and the scintillators. The technology is well established: a Shashlyk 
calorimeter used in E-865 at the AGS achieved S%/VE resolution, and ten 
years ago 75% of the PHENIX calorimeter was built with a similar design and 
shows similar performance. Subsequent studies done for the KOPIO experi- 
ment have demonstrated that a 3%/ \/E resolution is achievable. The designers 
of this calorimeter are now engaged in a design study of a Shashlyk calorimeter 
for STAR. 

In subsequent feasibility studies, a model of direct photon production has been 
constructed which agrees with both existing measurements at Pt > 2.0 GeV/c 
and a recent calculation below that value [28]. The model assumes Bjorken 
(boost invariant) expansion, and thermal spectra from three regions of tem- 
perature at different times in each rapidity slice. The ir^ spectrum is taken 
from PHENIX data. The rapidity distribution is known from the PHOBOS 
measurement, however, in the STAR range [r] = ±1.0) it is effectively con- 
stant. An analysis procedure for the calorimeter has been devised which takes 
into account both the effects of isolation cuts and of unavoidable unrecognized 
overlaps of background photons. 
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Fig. 40. Simulated correlation functions using 40 million events with full background, 
as described in the text. Top left: Qinv) top right: Qos\, bottom left: Qxyz- 

STAR can measure the direct photon yields and HBT parameters for px < 
600 MeV/c with about 40 milhon events. Fig. 40 shows the simulated results 
for three one- dimensional correlation functions: Qinv is the usual invariant 
four-momentum difference; Qosi is the four-momentum difference in the usual 
out-side- long Bertsch-Pratt system [177]; and Qxyz is the four- momentum dif- 
ference in the x, y, and z Cartesian system. The figure shows that, despite the 
very low A, quite good signals can be obtained. The individual out, side, and 
long radii are functions of the photon kx, described in detail in a forthcoming 
paper [176]. The number of events necessary to observe direct photon correla- 
tions for pt between 1.0 and 2.0 GeV/c is being studied, but it would clearly 
involve much larger data samples. An interesting possibility is to use the high 
Et trigger to access the higher pt region since it measures the system size at 
a different time [175]. 



4-2.2 Spectral Distributions of the ai Meson 

An experimental measurement of the iso- vector-axial vector (ai) spectral func- 
tion is desirable because of its important role in the search for chiral symme- 
try restoration. It is related to measurements of thermal low-mass dileptons, 
dominated by the p (iso- vector-vector correlator), the chiral partner of the 
Oi, see Sec. 2.7. It has been suggested in [129] to make such a measurement 
via the 7r^7 invariant-mass spectra. Calculations of photons from ai decays 
have concentrated on the inclusive photon yield. It is very difficult to isolate 
the ai contribution and determine the spectral shape of the ai from inclusive 
measurements alone. This is the main motivation for measuring the associ- 
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Fig. 41. The spectral shape of the initiany-produced ai — > 77r (blue line) compared 
to the spectral shape of the observed oi — > jtt which have passed through the 
medium without being absorbed (green hne), normahzed to the number of produced 
ai. 

ated charged pion even though they experience strong final-state interactions, 
including absorption. Therefore, the in-medium spectral information will pri- 
marily pertain to the more dilute stages of a heavy-ion reaction. However, the 
Oi may undergo significant modifications even at moderate temperatures and 
densities [178]. Such effects have possibly been seen in the p tt+tt^ [179] 
decay, which is "penalized" by two pion absorption factors. Important infor- 
mation about chiral symmetry restoration may thus be obtained together with 
the dilepton measurements of the p meson. 

We report on studies [180] of the ai —>■ 777 channel using transport model 
simulations, i.e., minimum bias Au+Au UrQMD events [181] at ^/s^ = 200 
GeV. In UrQMD, the ai was introduced with B{ai —>■ 777) = 0.1, B{ai —>■ 
pn) = 0.9, Ma, = 1230 MeV, and = 400 MeV. Only 5% of the produced 
ai's are not absorbed in the medium. Of these unabsorbed ai's, 80% decay 
via 77r and 20% via pre. The spectral shape of ai •yn that are not absorbed 
(green line in Fig. 41) is significantly different from the spectral shape of the 
produced ai's (blue line in Fig. 41). The peak is shifted toward lower masses 
by about 200 MeV. Thus, the extraction of medium effects seems feasible but 
requires a good understanding of final-state absorption effects. 

At RHIC-II, the ai — > 777 decay can be measured using the STAR full cov- 
erage TOF. In STAR, the efficiency of measuring a conversion 7 is 5% [182]. 
Therefore, 55M minimum bias Au-|-Au collisions are necessary to measure a 
3a signal (using HIJING to estimate the background). If a Shashlyk calorime- 
ter measuring low-momentum 7's could be used in STAR, the ai efficiency 
could improve by a factor of 10. 
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Fig. 42. Left: Measured vr*^ Raa as a function of for the — 10% most central 
Cu+Cu collisions at yJSNN = 22.4, 62.4, 200 GeV in comparison to a jet quenching 
calculation. The error bars in this figure represent the quadratic sum of the sta- 
tistical uncertainties and the point-to-point uncorrelated and correlated systematic 
uncertainties. The boxes around unity indicate uncertainties related to {Ncow) and 
absolute normalization. The bands for the theory calculation correspond to the as- 
sumed range of the initial gluon density dN^ /dy. The thin solid line is a calculation 
without parton energy-loss for central Cu+Cu at ^/snn = 22.4 GeV. Right: The 
average Raa in the interval 2.5 < pT < 3.5 GeV/c as a function of centrality for 
Cu+Cu collisions at ^/snn = 22.4, 62.4, and 200 GeV. The shaded bands repre- 
sent jet quenching calculations at three discrete centralities (A'^part ~ 10,50,100). 
The boxes around unity represent the normalization and (A'^coii) uncertainties for a 
typical A'coU uncertainty of 12 %. 

4-3 High Statistics and Energy /Species Scan 



The RHIC-II project calls for a ten-fold increase in the luminosity currently 
provided by RHIC, essential for full exploration of rare probes and very high-p^ 
signals that are currently statistics-limited. Examples include detailed heavy 
quarkonium spectroscopy, more precise values and pr-evolution of photon and 
hadron Raa at very high p^, as well as photon flow and HBT at high p^p. 
Photon-tagged jets, particularly if measured relative to the reaction plane, 
would give further insight to energy loss mechanisms along with the com- 
plementary jet-dilepton correlations. At the current luminosity, it would take 
several years to collect sufficient statistics for these measurements with just 
one species and one energy setting: clearly inadequate for exploring some of 
the key questions about the sQGP in depth and in a timely fashion. 

There are, however, additional physics opportunities. The above measurements 
investigate the properties of the sQGP above T^, certainly a necessary task. 
But, as emphasized in the Introduction, while the existence of a new state of 
matter has been reasonably well established at RHIC, key features of the phase 
transition itself, including variations across the T—^b phase diagram, have not 
yet been mapped out. The only way to do this is by a detailed energy/species 
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scan, only possible over a reasonable time scale with RHIC-II luminosities. 
The prospects are well illustrated by Fig. 42 which shows 7r° Raa in Cu+Cu 
collisions for 22.4GeV (close to the SPS range) to the highest RHIC energy, all 
measured by the same experiment (PHENIX). The transition from Raa > 1 
to strong suppression occurs in this energy range. Clearly this doesn't exclude 
the possibility of some suppression already at SPS energies, but the (Cronin) 
enhancement is still dominant - and essentially constant - at all centralities, 
whereas at 62.4 GeV the suppression increases with centrality. 

At current RHIC luminosities, even with the accelerator running for ~ 22 — 28 
weeks per year, the practical limit is at most two species or three energies per 
run period. The dominant part of the run is spent on actual collisions and 
data taking rather than changing settings. Since low px observables (bulk 
observables, thermal photons and dileptons, and hadron ratios) explore the 
phase transition, long runs are not required to obtain reasonable statistics 
for a given energy /species setting. Indeed, with the luminosity upgrade, suffi- 
cient data could be collected in a few weeks at most settings, even after the 
luminosity loss with 7"^ is accounted for. 

The baryon chemical pontential, hb, can be varied by changing -^/s^. RHIC 
was designed to vary ^/s^ and has successfully done so in several runs. The 
RHIC experiments have collected heavy-ion data at y/s^^ = 22, 62, 130 and 
200 GeV, covering several points in the y^s^ interval between the SPS and 
the maximum RHIC energy. Recent studies have shown that RHIC could run 
collisions as low as sqrtSj^j^ — 5 GeV — the AGS range — to compare to previ- 
ous AGS and SPS experiments, as well as fill in gaps in earlier measurements. 
While this may at first appear to be a repetition of earlier work, we want to 
re-emphasize the value of doing measurements over the entire energy range 
with the very same apparatus. Since most systematic errors are the same, the 
evolution of physics quantities characterizing the transition can be traced to 
higher precision and with more reliable systematics. 

Lower energies, however, come at the price of an approximately quadratic lu- 
minosity decrease with ^/s^. Therefore, as indicated at the end of Sec. 4.1.1, 
the RHIC-II luminosity upgrade is critical for a meaningful energy-scan pro- 
gram for low-mass dileptons. For example, the RHIC-II luminosity at y/sj^ = 
30 GeV will be the same as the current maximum-energy luminosity. Without 
upgrade, a factor of 40 or more luminosity loss from y^s^ = 200 GeV to 
^/s^ = 30 GeV would make an energy scan program essentially useless for 
dileptons. Based on discussions in Sees. 2.6, 2.7.2 and 3.1, a minimal ~ 20% 
accuracy with a ~ 2% mass resolution is needed to develop sufficient discrim- 
inating power between even rather distinct in-medium modification scenarios. 
Obviously, future theoretical refinements will further raise these requirements. 
Direct photon fiow measurements up to pt ~ 8 — 9 GeV are needed to inves- 
tigate contributions from parton fragmentation, jet-thermal interactions and 
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Bremsstrahlung off quarlcs. Current luminosities limit the measurements to 
~ 5 — 6 GeV, see Sec. 3.5. A tenfold luminosity increase would extend the 
PT-range by the necessary 2-3 GeV. Interferometry of pre-hadronic photons 
(emitted before the hadronic phase), pr > 2 GeV [175], with the precision de- 
duced in Sec. 4.2.1 requires about 7 x 10^ events which would, with current lu- 
minosities, allow measurements with only one system and/or energy per year. 
Also, many of the intermediate invariant mass and transverse momentum, 
M,pt ~ 2 GeV, signals, where sQGP radiation dominates (see Sees. 2.4.1 
and 2.5.3), can be triggered, making high luminosity extremely valuable in 
this context as well. Photon-based temperature measurements, described in 
Sec. 2.4.2, depend on detailed energy and species scans, including sufficient pp 
reference data at the same energies as A A collisions instead of relying on ear- 



lier results from other accelerators I I . Once again, the benefits of measuring 



both spectra in the same detector are obvious: in many cases the systematic 
errors decrease by 40 — 50%. 



Colliding different ion species allows variation of the collision geometry at 
fixed participant number, A^part- So far at RHIC several bulk observables have 
been found to scale with Apart irrespective of collision geometry. Recently it 
has been found that even more subtle quantities such as the azimuthally- 
integrated nuclear modification factor, Raa, is similar in central Cu+Cu colli- 
sions and mid-peripheral Au+Au collisions with the same Apart- However, once 
azimuthal distributions are studied, significant differences emerge since the two 
collision geometries are not alike (spherical Cu-|-Cu vs ellipsoidal Au-fAu) and 
these different geometries are crucial for understanding processes such as in- 
medium energy loss. Also, for sufficiently light ions, thermalization, apparent 
in the hydrodynamic behavior of the Au+Au and even Cu+Cu data, can no 
longer be achieved. Pinpointing the thermalization boundary could be cru- 
cial for determining the behavior of probes near the phase boundary and the 
conditions needed for the transition to occur. 



We reiterate that higher luminosities can be used to increase the statistics of 
a particular data set, gaining access to rare probes, and also to increase the 
number of ion species (symmetric or asymmetric collision systems) and/or the 
energy settings by subdividing the RHIC runs since many important signals 
do not require very high statistics. The accelerator is the most flexible ever 
built and can run multiple system/energy combinations per year if the physics 
justifies it. In fact, RHIC produced the very first collisions at y/s^ = 9GeV 
in a 12 hour test in June 2007. 
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Sec. 3.7 
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Sec. 2.4.1 
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V2 
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n, 
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Sec. 4.3 
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pQCD, scale 
jet energy scale 
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Sec. 4.3 
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Sec. 4.2.1 
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Sec. 3.5 
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ai mass 
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absorption? 


Sec. 4.2.2 
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rapidity 


7- jet 

Xc 


CGC 


forward 
calorimetry 


common 7, tt^ 
suppression? 


Sec. 4.1.3 



Table 2 

Bottom-up view of electromagnetic probes in heavy-ion collisions: from observables 
to the pertinent theoretical issues addrssing the underlying physics questions. 

4-4 From Observables to Pertinent Theoretical Issues 

In Table 2 we summarize how electromagnetic probes can influence the theo- 
retical understanding of relativistic heavy-ion collisions. In the first column we 
list specific measurements, in the second the signal extracted either directly or 
with minimal model assumptions, and the physics question addressed in the 
third column. The fourth column describes the most important experimental 
issues and the fifth indicates open or ambiguous theory questions that require 
significant further investigation. Finally, the sixth column points either to a 
Section of this review or to a reference where the issues are discussed in more 
detail. 

See Sec. 3.4 for an example of the difference this might make. 
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While this tabic has a level of subjectivity, it reflects our current thinking on 
the most interesting issues related to EM probes at RfflC-II. We do not list 
topics that are primarily or exclusively the domain of the LHC or the future 
Electron-Ion CoUider (EIC). 



5 SUMMARY AND CONCLUDING REMARKS 

We now summarize our opinion of the most promising developments in con- 
nection with EM probes at future RHIC experiments. 

Low-mass (axial-) vector-meson spectroscopy close to the chiral transition 
Precision dilepton data will provide detailed information on medium modifi- 
cations of p, cu and (f) mesons, and thus illuminate the question of hadronic 
mass (de-)generation. "Quality control" of theoretical models via independent 
constraints from symmetries, QCD sum rules and phenomenology is essential 
to limit the scope of viable axial/vector (A/V) spectral functions. An exci- 
tation function measurement will discriminate temperature and net-baryon 
density effects to systematically map out in-medium effects across a signifi- 
cant part of the QCD phase diagram. Information on chiral symmetry should 
be inferred from a novel measurement of the axial vector spectral function 
via 7r^7 invariant-mass spectra as well as through a well-defined combination 
of effective chiral hadronic models, chiral sum rules and finite-temperature 
QCD lattice computations of chiral order parameters such as the pion decay 
constant and four-quark condensates. 

The required theoretical tools are largely in place: chiral effective models for 
realistic axial-/vector spectral functions at finite temperature can be used to 
calculate the temperature dependence of order parameters (moments of "V — 
A" spectral functions). Unquenched IQCD evaluations of the latter should be 
pursued with high priority. A convincing deconvolution of the vector correlator 
from the measured spectra will require reliable space-time descriptions of AA 
collisions, expected to emerge from progress in hydro dynamical simulations 
and complementary transport models. 

Footprints of chiral restoration are furthermore expected at dilepton masses 
around the ai mass, 1 — 1.5 GeV, due to "chiral V — A mixing". To detect a 
continuum enhancement of less than a factor of two requires accurate charm 
and background determinations providing a signal with no more than 20% to- 
tal error, achievable with the planned vertex-detector upgrades together with 
sufficient statistics. In addition, dilepton elliptic flow and transverse momen- 
tum spectra will provide additional information on the space-time emission 
profile at different invariant masses. 
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The highest temperatures of the matter formed at RHIC 

We have identified three promising regimes where the QGP is expected to be 
the dominant source of electromagnetic radiation: intermediate-mass dilep- 
tons (1.5 < M < 3 GeV) and photons and low-mass dileptons at intermediate 
transverse momentum, (1 < ?t < 3 GeV). These measurements will be sup- 
plemented by 77 correlation analyses which yield complementary information 
on the temperature and fireball size. 

QGP resonances 

The only direct way to experimentally determine the existence of hadronic 
bound states/resonances in the sQGP is to search for a resonant dilepton 
signal. Due to the nature of QGP emission, the greatest sensitivity to perti- 
nent vector states is in the M = 2 GeV mass region, coinciding with current 
expectations from lattice QCD and effective models. RHIC provides optimal 
conditions for this search since the initial temperatures are close to the antic- 
ipated bound state dissolution temperatures. 

As is the case for signatures of chiral mixing at lower masses, an accurate 
determination of background and charm sources is mandatory to restrict the 
total experimental error to below the ~ 20% level. If the resonance states exist 
at masses below ~ 1.5 GeV, discriminating them from chiral-mixing effects 
will be more involved, increasing the demand for accuracy and systematic 
centrality and excitation-function studies. 

Detector and luminosity requirements 

To achieve the required background rejection and obtain high-precision charm 
measurements, PHENIX needs the Hadron Blind Detector and vertex detec- 
tor, while STAR needs the Time-of-Flight detector. Heavy Flavor Tracker and 
Data Acquisition System. Initial dilepton measurements can be made after 
these new detectors are in place. However, an energy scan with statistics com- 
parable to recent NA60 data requires significant luminosity upgrades. 

A summary is provided in Table 2. 

Based on the arguments given in this review, we believe that after the fu- 
ture RHIC detector and luminosity upgrades for heavy-ion collisions in the 
10 < ^/s^ < 200 GeV regime, electromagnetic probes will, in a combined 
experimental and theoretical effort, result in decisive and unique new perspec- 
tives of QGD matter at high (energy-) density. In particular, new insights into 
hadronic mass generation mechanisms related to chiral symmetry restoration, 
thermal radiation at unprecedented temperatures, microscopic properties of 
the sQGP and system size measurements during the early matter evolution 
will be gained. 
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